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ABSTRACT 

This report describes an interpretive system which transforms the 650 into 
a three-address, floating-decimal, general-purpose computer, primarily suited 
for scientific and engineering calculations. The system is complete in the sense 
that all mathematical, logical and input-output operations normally called for in 
such calculations can be performed within the system, i. e. , without reference to 
the basic operation codes of the 650. The guiding principles in designing the sys- 
tem have been ease of use, as defined in the introduction, high speed of arith- 
metic and frequently used logical operations and full accuracy and range for the 
elementary transcendental functions. 

The report serves a dual purpose. It presents the external characteristics 
of the interpretive system to the potential user by means of detailed explanations 
accompanied by illustrative examples, assuming no previous familiarity with 
internally programmed machines. It also describes the internal structure of the 
system to the professional designer of such systems, enabling him to modify it 
to suit his particular needs or to borrow ideas or building blocks he may find 
useful. 

The system is available in punched card form to anyone who requests it. 



-1- 



CONTENTS 



Note: The material of immediate concern to those who wish to learn how to 
program problems in the interpretive system is contained in sections II-X. Sec- 
tion I is devoted to general considerations and may be bypassed. Section XI deals 
with the internal structure of the system, primarily for the benefit of those inter- 
ested in the design of interpretive systems, but the discussion of possible mod- 
ifications in Sec. XI. 1. and the contents of Sec. XI. 2. and XI. 3. should be of wider 
interest and do not require familiarity with the basic language of the 650. 

The experienced programmer may absorb the essentials of the system by 
reading the definitions of the operations. Page references to them are given in 
the summary of operation codes. 



Page 



I. INTRODUCTION 4 

I. 1. General Design Considerations 4 

I. 2. Changes and Additions 6 

II. GENERAL INFORMATION 7 

II. 1. The 650 7 

II. 2. The Interpretive System: 

Storage; Data and Instruction Form 7 

III. MATHEMATICAL OPERATIONS 10 

III. 1. Arithmetic Operations 10 
III. 2. Special Functions 10 
III. 3. MOVE 000 11 

III. 4. An Example 12 

IV. LOGICAL OPERATIONS 13 

IV. 1. Transfer Operations 13 
IV. 2. Loop Operations 15 
IV. 3. Address Change Operations 18 

IV. 4. MOVE 23 

V. INPUT-OUTPUT OPERATIONS 25 

V. 1. Card Form 25 



CONTENTS (continued) 



Page 

V. 2. Punching 25 

V. 3. Loading 27 
V.4. Reading 29 

VI. PROGRAM TESTING 31 

VI. 1. Memory Print-Out 31 



VI. 2. Tracing 



32 

VI. 3. Console Testing 34 

VII. SUMMARY OF OPERATION CODES 36 

VIII. STOPS 37 

IX. EXECUTION TIMES 40 

X. SPECIAL, TOPICS IN PROGRAMMING 43 

X. 1. Subroutines; Translation 43 
X. 2. Unnormalized Input; Transition Between 

Floating- and Fixed-Decimal Form 44 

X. 3. Examples 44 

XI. INTERNAL STRUCTURE OF THE SYSTEM 50 

XI. 1. Detailed Design Considerations 50 
XI. 2. Related Systems 55 
XI. 3. Numerical Methods 55 
XI. 4. Control Panel Wiring for the Read-Punch Unit 57 
XI. 5. Control Panel "Wiring for the Tabulator 58 
XI. 6. System Loading 6l 
XI. 7. Programs 63 



•3- 



I INTRODUCTION 



I. 1. GENERAL DESIGN CONSIDERATIONS 



The use of most existing computing devices whose degree of automatic per- 
formance substantially exceeds that of a desk calculator entails certain problems 
not encountered in desk computing. To cope with these problems, one may incor- 
porate additional circuitry into the machine --this, indeed, appears to be the trend 
in recently announced commercially available machines- -or, alternatively, one 
may program, in terms of the basic language of the machine, a system or super- 
language in terms of which the general user will program his problems. The user 
may consider the machine and the super -language as one entity, -and no knowledge 
of the basic machine language is required of him. Before actual calculation, the 
programmer's instructions are translated by the machine into the basic language. 
If this translation or interpretation takes place each time an instruction is to be 
executed, rather than once for all at the beginning of a problem, the super-language 
is referred to as an interpretive language or system. Limitations in storage capac- 
ity may necessitate the choice of an interpretive system rather than a system of 
the once-for-all type in the case of most small or medium- sized computers. 

The designers of an interpretive system are faced with a very large number 
of decisions. To provide a basis of motivation for these decisions, it is convenient 
to list here, in somewhat arbitrary order, some of the above-mentioned problems 
which the present interpretive system proposes to solve. All of them may funda- 
mentally be measured in terms of total time spent by a programmer in learning to 
use the machine and in using it on a specific problem. In this sense, the "ease of 
use" referred to in the abstract above is implicitly defined by the list that follows. 
The price paid for the saving of programmer time is, of course, to be found in 
substantially reduced speed of operation. 

A. Scaling 

The storage medium- -paper- -normally used in desk computing places no 
practical restriction on the size of numbers or on the location of decimal points. 
In using a computer that automatically stores intermediate results in registers 
of fixed length and with the position of the decimal point fixed in advance, a great 
deal of time must in most cases be spent on estimating the range of all intermediate 
results to prevent errors due to overflow. The well-known way of avoiding this 
at the expense of a very substantial increase in the internal complexity of the arith- 
metic operations is to represent each non-zero number in floating -decimal form, 
i.e., as a signed quantity whose absolute value lies in a fixed range, accompanied 
by an exponent of 10 or decimal point indicator. 

B. Length and Complexity of the Program 

Floating decimal arithmetic and frequently needed special functions could 



-4- 



be incorporated into a program written in the basic machine language in the form 
of a set of subroutines reached by a two-way transfer of control, ("calling se- 
quence" or "basic linkage") and there are indeed problems for which this is the 
best choice. In many cases -particularly in the case of relatively short prob- 
lems where the results are needed quickly- -a further reduction of the program- 
ming effort is desirable. This may be achieved by combining under single oper- 
ation codes, in an interpretive language, groups of steps in the bas 1C language 
needed for performing frequently occurring tasks. For example, a single instruc- 
tion in which three locations are specified may be used for adding two arbitrarily 
located numbers and storing the result or a block of information of any length 
may be punched on cards as a result of a single instruction. In particular, the 
task of repeating a calculation a specified number of times, each time with appro- 
priate modifications, must be made easy, and the interpretive system described 
in this report goes as far as is believed possible in this direction by providing an 
order ("LOOP") with which simple cases of this task can be handled by a single 
instruction. 

To preserve the simplicity gained by introducing an interpretive system, the 
system must be made complete or self contained so that most problems can be 
conveniently programmed without reference to two different systems of °P eratl ™ 
codes, although, of course, leaving and re-entering the interpretive system should 
be made possible in order to provide the experienced programmer with complete 
flexibility. 

C. Restrictions 

In desk computing, one cannot fail to notice if the argument for which a func- 
tion value is to be found in a table falls outside the range covered by the table or 
if in transcribing a set of numbers from one area on a piece of paper to another, 
overlapping area, some of the numbers to be transcribed are erased before tran- 
scription. A machine will avoid or detect such blunders only if programmed to do 
so, and as much as possible of this programming should be included m the inter- 
pretive system. For example, it is desirable that all mathematical f^ c *™f «- 
cluded in the system be available for the full range of argument consistent w lt h their 
definition and with the limitations imposed by the machine itself and that they are 
computed to the full accuracy of the number system used. Error stops «*«£* 
violations of unavoidable restrictions should be included to the fullest extent that 
space limitations permit. 

D. Program Testing 

The usefulness of a general -purpose computer or interpretive system depends 
decisively on the methods provided for testing ("debugging") new programs-for 
definitions and details, see Sec. VI. In the case of an interpretive system whose 
operating epeed'is only one order of magnitude above the speed of card pmch^g, 
testing by means of a tracing routine included in the system compares favorably to 
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console testing, at least in the case of programmers whose familiarity with the 
machme is limited. Either of these methods is thought of as a tool normally used 
only when memory print-outs have been found insufficient. To facilitate testing 
by any of the methods mentioned and to keep programming as concrete as possible, 
the system described in this report assumes that the actual machine location of 
each instruction and stored number is assigned by the programmer. The system 
may, of course, be used in conjunction with regional or symbolic assembly programs 



I. 2. CHANGES AND ADDITIONS 



«. r ^ n ?l° US mm ° r chan S es su §g est themselves when the system is viewed in 
the light of the experience gained in designing it; some of them are discussed in 
detail at the end of the report. Corrections of errors not yet revealed must be ex- 
pected. External changes and additions will undoubtedly be proposed after a period 
of use. The present system should, therefore, be considered primarily as a first 

Ind'ZJt eaCh K US I rmay C ° nSider chan S in S to bett <~ «uit his needs. Comments 
and suggestions on both internal and external aspects of the system will be greatly 
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II GENERAL INFORMATION 



II. 1. THE 650 

The IBM 650 is an electronic computer whose basic storage consists of 
a magnetic drum capable of holding 2000 words (numbers) of ten decimal digits 
and sign. The machine is internally programmed, i. e. , the program of instruc- 
tions which the machine is expected to follow is kept on the drum, and the ma- 
chine automatically reads one instruction at a time from the drum, executes it, 
reads another instruction, and so on. Initially, the program is loaded onto the 
drum from punched cards but each instruction is loaded only once, although the 
machine may be expected to execute it many times in the course of a problem. 
Special orders are inserted into the program to cause repetition of prescribed 
sections the desired number of times. In many cases, the instructions exec- 
uted by the machine are changed or modified under control of the program between 
successive executions. This ability to modify its own instructions is one of the 
distinguishing characteristics of an internally programmed machine. 

In the basic model of the 650, all answers are punched by the machine into 
cards, which may be printed on a separate tabulator. 

The 650 is a general-purpose, fixed-decimal machine and any programmer 
may, with the aid of a detailed manual published by IBM, learn to use it as such. 
There are many large problems and many problems of a data processing nature 
for which fixed-decimal operation is definitely indicated, and the programmer is 
asked to give serious consideration to this alternative for all but the very smallest 
problems, since the gain in machine time over floating -decimal operation (ex- 
plained below) may be as high as 10:1. The machine -language programmer may 
relieve himself of many tasks by using the interpretive system for loading, punch- 
ing, calculation of special functions, etc. , provided 1000 storage locations suffice 
for his problem. (See TR OUT, Sec. IV. 1) 

This report describes a system which enables the programmer to use the 
650 as a floating -decimal machine, without being familiar with the fixed-decimal 
mode of operation. Beginning with the next section, all statements will concern 
the system rather than the 650 itself, but it should be borne in mind that if any- 
thing in the system appears restrictive from the viewpoint of a particular appli- 
cation, —storage capacity, speed, card form, word length, etc. --total or partial 
use of basic 650 coding may be the answer. 

II. 2. THE INTERPRETIVE SYSTEM: STORAGE; DATA AND INSTRUCTION FORM 

When the interpretive system is in use 999 ten-digit storage locations, 
numbered 001-999, are unrestrictedly available to the programmer for storing 
instructions and data. The location 000 has a special use ("previous result") 
which will be discussed below. 
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Throughout the system, numbers upon which mathematical operations are 
performed are stored and used in so-called (normalized) floating -decimal form, 
which will be defined as follows: The number zero is written as ten zeros with a 
plus sign ("machine zero"). Any number A other than zero is expressed as 



A= ± A, 



10 a l 



where 1^ Aj<. 10 and -SOj^a^ 49. In the machine, A is written as the pair 
±(Aj, a), where a = aj+ 50 and Ai is an eight-digit number with seven decimal 
places. The "machine exponent" a is a two-digit (positive) integer located at the 
right end of the number. Non-zero numbers A not in the range 10" 50 ^ | Al<10 50 
cannot be correctly used in the machine, and some of the mathematical operations 
will give an error stop if the result would fall outside this range (see STOPS). Num- 
bers loaded into the machine must also be in the form prescribed above, unless 
special precautions are taken (see Sec. X. 2). 



The system instructions are signed ten-digit numbers of the following 



+ 


°1 


A or 2 


B 


C 



form: 



Here, X is a one-digit operation code and B and C are three-digit addresses. The 
three-digit quantity "A or 2 " is interpreted as an address A if X ^ and as an 
operation code 2 if 1 = 0. The sign of the instruction is used in connection with 
the LOOP order (see LOOP OPERATIONS). The only difference between the mu- 
tually exclusive L and 2 operations is that all operations which require three 
addresses have been designated 0^, all others, O2. 

An example will illustrate how the addresses are used in a program: To add 
the number stored in register 200 to the number stored in register 201 and store 
the result in 500, the operation code L = 1 (ADD) is used and the instruction reads: 
1 200 201 500. To take the square root of the number in 200 and store the result 
in 500, the 2 -operation 2 = 300 (SQRT) is used: 300 200 500. As will be shown 
in later sections, it is also possible to call out an instruction stored in memory 
and operate upon it, e.g. , increase one of the addresses in it. In storage, no dis- 
tinction is made between instructions and data so that the programmer is free to 
use any memory location for storing an instruction or a number as he sees fit. 

To facilitate explanations, the following notation will be used: The ten-digit 
quantity jvhose storage location has the address A_will be denoted by A; anal- 
ogously, B will denote the contents of location B. C denotes the result of a calcula- 
tion about to be stored in location C. 

In addition to being stored at C, the result C of any mathematical operation 
(i. e. , arithmetic operations or special functions) and of MOVE 000 and CONS 
(see READING) is automatically stored in the special location 000. If this result 
is needed on the next step, (or, more generally, before it has been replaced by 
the result of a subsequent mathematical operation or MOVE 000 or CONS) calling 
it out by using 000 as an A-address will reduce the execution time in the case of 
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the arithmetic operations. Also, time will be saved in any mathematical operation 
by using 000 as a C-address when C will be needed only on the next step. Execu- 
tion times are discussed in detail in a later section, but it should be emphasized 
that timing considerations only affect the running time of a problem, never the 
correctness of results. All locations are accessible at any time. 

Special addresses for obtaining frequently needed numbers, such as, and 
1, are not provided by the system. The programmer should load such numbers 
into locations of his own choosing. 



Ill MATHEMATICAL OPERATIONS 



III. 1. ARITHMETIC OPERATIONS 

The operations will be introduced in an order chosen to' facilitate the learn- 
ing process. Later, a concise summary of all operation codes will be given. 
Alphabetic operation codes are listed in addition to the numerical ones merely to 
facilitate programming; they are not introduced into the 650 and need not be used 
at all. 

The result of each arithmetic operation is rounded. If the result is zero, 
a machine zero is given, i. e. , the machine exponent will be 00. An error stop 
occurs if the result of a multiplication or division would fall outside the range of 
the floating -decimal number representation; another error stop detects attempts 
to divide by zero (see STOPS). 



A list of the arithmetic operations follows: 



Numer. 
code 

0, = 1 



Alpha, 
code 

ADD 



Function 



0, = 2 



SUB 
MPY 
DIV 

NGMPY 
III. 2. SPECIAL FUNCTIONS 



°1= 3 



0, = 4 



'1 



°1= 5 



Add jin floating -decimal form) the num- 
ber A stored at A to the number B 
stored at B, store the result C at C 
and 000. Abbreviated: 

A+ B = C 

Subtract: A - B = C 

Multiply: A . B = C 

Divide: A / B = C 

Multiply negatively: -A . B = C~ 



The system is intended to give eight-digit accuracy (i. e. , an error less than 
1 in the eighth digit) in computing the special functions included whenever the input 
makes this accuracy possible. For trigonometric functions of an argument exceed- 
ing one revolution and for logarithms of numbers near 1, loss of accuracy follows 
from the mathematical properties of the respective functions and stops (which may 
be bypassed by the setting of a console switch) are provided when this loss exceeds 
two digits. For small values of the argument, where an eight-digit, fixed-decimal 
representation of the sine or arc tangent would contain leading zeros, the floating- 
decimal representation would normally introduce meaningless digits at the right 
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end. To reduce this nuisance to a tolerable level and also make possible trig- 
onometric calculations with extremely small arguments, the formulas sin x = x 
and arc tan x = x are used for |x\< . 0025 and |x|< . 001, respectively. Those 
interested will find the methods of computing the special functions described in 
Section XI. 3. 

Aside from the limitations imposed by the above mentioned inherent loss 
of accuracy and by the floating -decimal representation of the result, no restric- 
tions apply to the natural range of the argument for the special functions. Error 
stops will prevent attempts to take the square root of a negative number or the 
logarithm of a non -positive number. 

The special functions (or, more precisely, elementary transcendental 
functions) are: 

Numer. Alpha. Function 

2 = 300 SQRT 71" = C 

2 = 301 EXP E e B = C 

2 = 302 LOG E log e B = C 

= 303 SIN R sin B = C, B in radians 

2 = 304 COS R cos B = C, B in radians 

= 305 ART R arc tan B = C, C in radians, ICKTr/2 

= 350 ABS |Bl = C 

2 = 351 EXP 10 10 B = C 

2 = 352 LOG 10 log 10 B = C 

= 353 SIN D sin B = C, B in degrees 

0-, = 354 COS D cos B = C, B in degrees 

= 355 ART D arc tan B = C, C in degrees, fc|< 90 

Subdivisions of a degree are expressed decimally, not in minutes and sec- 
onds. 

III. 3. MOVE 000 

In many cases (particularly in connection with the use of subroutines) it may 
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be convenient to be able to call out a number B from B and deposit it in C, as 
well as in 000, without the time- consuming use of a floating -decimal arithmetic 
operation. This is accomplished by the logical operation 0=9 ("MOVE") with 
A 5s 000. The normal use of MOVE with A £ 000 is described in Sec. IV. 4. 

III. 4. AN EXAMPLE 

For the benefit of anyone with no previous computer experience, a simple 
example illustrating the use of the mathematical operations will be inserted here. 
Suppose that, as a part of a program which is assumed to be already on the drum, 
it is desired to evaluate the function 



f(x) = 



sin x 



/ 



1+ e 



-x3 



Here, x in radians is assumed to be in storage register 500 and the constant 1 in 
600. The quantity e~ x3 is to be stored in 501, and f(x) in 502. A program might 
look as follows: 

C Comments 

000 x 2 

000 -x 3 

3 
501 e~ x , store in 501 

000. l+e- x3 



Alpha. 


°1 


MPY 


3 


NGMPY 


5 


EXP E 





ADD 


1 


SQRT 





SIN R 





DIV 


4 



or °2 


B 


500 


500 


000 


500 


301 


000 


000 


600 


300 


000 


303 


500 


000 


400 



400 ^1 + e" x3 , store 

temporarily 

00£ sinx 

502 f(x), store in 502 

The extensive use of the "previous result" address, 000, is worth noting. 
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IV LOGICAL OPERATIONS 



IV. 1. TRANSFER OPERATIONS 

Suppose the machine has been instructed (see LOADING) to begin a program 
by executing the instruction stored in, say, location 101. When this execution is 
completed, the machine will automatically execute instruction 102, then 103, etc., 
until told by the program to do otherwise. Operations whose primary function is 
to influence either the order in which instructions are executed by the machine or 
the selection of stored data upon which the instructions make the machine operate 
will be called logical operations. A simple example of such an operation is O2 = 
203, "Transfer Control". If in the sequence 101, 102, 103 above, instruction 103 
should read "TR 203 000 080", the next instruction executed by the machine would 
be 080 instead of 104. This may be expressed by saying that "control was trans- 
ferred to 080". The B-address was ignored in this case. The transfer of control 
may be made to depend on the result of calculations (mathematical or logical) in 
which case a "conditional transfer" is said to occur. Logical operations—con- 
ditional or unconditional- -are needed whenever several blocks of instructions, 
located on various parts of the drum, are to be tied together to form a program, 
whenever it is desired to repeat a calculation several times, etc. 

For simplicity in grouping, the following list of transfer operations includes 
two (UNC STOP and NOOP) whose transfer function is of a degenerate nature. In 
a first reading, it may be advantageous to omit the TR SUBR and TR OUT operations. 



Numer. 



2 = 000 



Alpha. 
UNC STOP 



Function 

Stop unconditionally. The machine 
stops regardless of the setting of 
console switches (see CONSOLE) 
and displays 9999 on the address 
lights and B on the display lights. 
This operation should be used only 
where it is intended to discontinue 
the execution of the program, since 
a continuation of the program cannot 
be effected by a simple depression 
of the PROGRAM START key (see 
STOPS). The C-address is ignored 
but should be filled, e. g. , with zeros 
(see LOADING). 



°2 = 



200 



COND STOP 



Stop conditionally and transfer. The 
machine stops if the PROGRAMMED 
STOP switch on the console is in 
the STOP position. The number 
1120 is displayed on the address 
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Numer. 



Alpha. 



Function 



lights and B on the display lights. 
When the PROGRAM START key 
is depressed, control is trans- 
ferred to C. If the PROGRAMMED 
STOP switch is in the RUN position, 
control is transferred to C with- 
out stopping. 

Caution: If the PROGRAMMED 
STOP switch is on RUN, the stops 
for loss of accuracy in sine, co- 
sine and logarithm and the stop in 
the CONS operation will not occur. 

This operation may be used for stopping at check points in the early running 
stages of a problem, with the option of avoiding the stops during later runs. 



0-> = 201 



TR SGN 



2 ss 202 



TR EXP 



Transfer on sign. Control is trans- 
ferred to C if the result of the last 
mathematical operation or MOVE 
000 or CONS is negative, to B if 
it is non-negative, (i. e. , zero is 
regarded as having a plus sign). 

Transfer on exponent. The exponent, 
c, of the result of the last mathe- 
matical operation or MOVE 000 or 
CONS is compared to B (the leading 
digit of B should be 0). Control is 
transferred to C if c > B. If c< B, 
control proceeds to the next instruc- 
tion. 



This operation is particularly suited for the summation of series where terms 
are to be added until they have a prescribed number (50 - B) of leading zeros. For 
example, in order to return to instruction 080 only as long as the absolute value of 
the previous result is . 0001 or greater, one would write "TR EXP 202 046 080". 
This saves a time-consuming floating subtraction preceding the test. The TR EXP 
operation is also intended to take the place of the TR ZERO operation found in most 
systems. Due to the accumulation of small errors during a calculation, it is unwise 
in most cases to expect a result to be exactly zero to eight figures; here a TR EXP 
with a suitably chosen B may prevent a never-ending repetition of a part of a program. 



2 ss 203 



TR 



Transfer. Control is transferred to 
C, i. e. , the next instruction exe- 
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Numer. 



Alpha. 



2 = 204 



TR SUBR 



2 - 205 



TR OUT 



Function 

cuted will be the one stored at C. 
The B address is ignored but 
should be filled, e. g. , with zeros. 

Transfer to subroutine. The C- 
address of the instruction located 
at C is set equal to B, whereupon 
control is transferred to C. The 
sign of the instruction at C is made 
positive. (For an elucidation and 
applications see SUBROUTINES. ) 

Transfer out. Control is transferred 
to C and the instruction stored there 
is executed in the basic language of 
the machine (i. e. , outside the inter- 
pretive system). When an instruc- 
tion address 1095 is given in the 
basic language, control is returned 
to the interpretive system beginning 
at the instruction following the TR 
OUT. The B-address of the TR OUT 
instruction is ignored but should be • 
filled, e.g. , with zeros. The pro- 
grammer in basic language must be 
careful not to use locations above 999, 
which are occupied by the interpretive 
system. 

No operation. Control proceeds to 
the next instruction. The B- and 
C-addresses are ignored but should 
be filled, e.g., with zeros. 

This operation is likely to occur in connection with tracing (see TRACING, 
particularly the ST TR ERAS operation) or when a superfluous instruction has 
been deleted from a program. 

IV. 2. LOOP OPERATIONS 

A highly repetitive character is required of any problem to be economically 
handled on an automatic computer. In certain instances, such as Newton's iteration 
procedure for the solution of equations, a repetitive process or "loop" is con- 
veniently programmed, merely using conditional transfer operations. In many cases, 
however, some of the instructions to be repeated must be slightly modified in a 
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systematic way before each new repetition. For example, in the evaluation of a 
linear expression £ a i x i with the aj and the x^ stored in blocks of consecutive 

i= 1 
locations, the addresses of a i and x^ must be increased by 1 each time a new term 
is to be computed. To facilitate programming of this kind, the system provides 
two methods of so-called address modification. The simpler — but less general-- 
of these methods employs a special counter called the loop box, which is stored 
in a location normally inaccessible to the programmer. If an instruction carries 
a minus sign, the current contents of the loop box will be added to the instruction 
(in fixed-decimal arithmetic and without regard to the sign) before it is executed. 
If, for example, the instruction - 1 531 600 901 is given and the loop box contains 
+ 009 000 009, the instruction actually executed by the machine would read 1 540 
600 910. The original instruction remains unchanged in its storage location. At 
the end of a calculation, an 2 instruction called LOOP enables the programmer 
to increase the contents of the loop box by 1 in one or several address positions and 
to transfer control back to the beginning of the calculation. Hence, the calculation 
may be carried out repeatedly, each time with different addresses used in the exe- 
cution of instructions with minus signs. A test provision included in the LOOP order 
stops the repetition after a specified number of executions and resets the loop box 
to zero for future use. An example will be given after the following list of LOOP 
operations. 



Numer. 



2 = 100 



Alpha. 
LOOP A 



Function 

Loop on A. The contents of the 
loop box are increased by 001 
000 000. In other words, the A- 
segment of the loop box is in- 
creased by 1. After the increase, 
the A-segment of the loop box is 
compared to the B -address of the 
LOOP instruction. If the A-seg- 
ment is less than B, control is 
transferred to C. If the A-seg- 
ment is equal to B, (or greater, 
which will never be the case in 
normal use) the loop box is reset 
to zero and control proceeds to 
the next instruction. 



0, = 010 



LOOP B 



Loop on B. Analogous to LOOP A, 
with the B- segment of the loop 
box now being increased and com- 
pared to the B -address of the LOOP 
instruction. 



0, = 001 



LOOP C 



Loop on C. Analogous to LOOP 
A, with the C- segment of the 
loop box being increased and 
compared to B. 



-16- 



Numer. 



2 = 110 



Alpha. 
LOOP AB 



Function 

Loop on A and B. Analogous to 
LOOP A. The A- and B -segments 
of the loop box are increased by 
1 and the A-segment is compared 
to B. 



o 2 = 101 



0? =r Oil 



°2 = 



111 



LOOP AC 
LOOP BC 

LOOP ABC 



Loop on A and C. Analogous. 

Loop on B and C. Analogous. The 
B -segment is used for the com- 
parison. 

Loop on A, B and C. Analogous. 
The A-segment is used for the 
comparison. 



To illustrate the use of a LOOP order, consider the evaluation of the linear 
expression L(x) = £? a- x-, where the a. and the x- are stored in memory. In 

ih x 

choosing storage locations for numbers, it is wise to plan in advance how they are 
to be used in the program. In this case, since the a i and the x. are to be reached 
using the LOOP operation, it is advantageous to store them in blocks of consecutive 
locations, say the a i in 800 + i and the x i in 900 + i, (i = 1, 2, .... 20). Suppose 
L(x) is to be stored in 700. For simplicity, assume that register 700 contains zero 
at the beginning of the calculation and that the loop box has been reset. The entire 
program for this calculation might be written as follows: 



Instr. No. 

101 
102 
103 
104 



Alpha. 



Sign. 







1 



A or 0, 



MPY — 3 

ADD + 1 

LOOP AB + 

Next instruction in the problem. 



801 


901 


000 


000 


700 


700 


110 


020 


101 



Note that the B-address of the LOOP order simply indicates the number of times 
the arithmetic calculation is to be performed, including the first time when the 
addresses are actually unmodified (modified by adding zero). The practice of 
starting the instruction numbering at, e.g. , 101, rather than 001 facilitates later 
additions to the beginning of a program. 

The loop box is automatically reset at the beginning of a new problem (see 
LOADING), and whenever a transfer out of a loop is effected by a loop order (as 
stated in the definitions above). Hence, the resetting of the loop box need not 
concern the programmer under normal conditions. If the need for resetting the 
loop box should arise, however, this is easily done by giving, e.g. , the order 
LOOP A with the B-address 000. According to the definition of LOOP A, this will 
cause control to proceed to the next instruction with a resetting of the loop box. 
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The C-address is irrelevant in this case. This situation would arise if control 
were transferred out of a loop in the middle of it by one of the conditional transfer 
operations. 

It is worth observing that a LOOP operation may be advantageously used in 
some cases where address modification is not involved, simply to repeat a sequence 
of steps a prescribed number of times, e. g. , each time adding a fixed increment to 
a parameter. In such a case, any one of the loop orders could be chosen, (see 
EXECUTION TIMES, however) and no negative instructions would occur. 

The advantages of the loop-box method are its simplicity and high speed and 
the fact that the original instructions remain unchanged in memory. It is limited 
by the fact that there is only one loop box and hence, all instructions to be modified 
are modified in the same way. To handle situations more complicated than this, 
the system provides a set of operations described in the next section. 

IV. 3. ADDRESS CHANGE OPERATIONS 

Many problems can be completely programmed without the use of address 
change operations, and for someone approaching the field of internal programming 
for the first time, it might be advantageous to ignore this section until the need 
for more general logical operations arises. 

The functions of the address change operations are: (a) To increase or 
decrease a designated address of an instruction in storage by any given amount; 
(b) To set such an address to a given value (without reference to its previous value); 
and (c) To transfer control as a result of comparing such an address to a given num- 
ber. 

There are nine O^-operations among the address change operations. In each 
of these, the B-address gives the location of the instruction (B) to be changed and 
the C-address is the amount of change. For example, suppose the instruction 
600 750 005 (using the operation 2 = 600,. ADD A) is given and suppose location 
750 contains the instruction I 320 400 000. Then the A-address, 320, of this in- 
struction will be increased by 005 and the resulting instruction 1 325 400 000 stored 
back in 750. Similarly, if 050 750 333 were given, (using 2 = 050, SET B) the 
instruction in 750 would be changed to read 1 320 333 000. In brief: 

Numer. Alpha. Function 

2 = 500 SET A Set the A-address. _ The A-address 

of the instruction (B) specified 
by B is set equal to C. 

°2 = 05 ° SET B Set the B-address. _ The B-address 

of the instruction (B) specified by 
B is set equal to C. 

°2 = 005 SET C Set the C-address._ The C-address 

of the instruction (B) specified by 
B is set equal to C. 
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Numer. 



2 = 600 



Alpha. 
ADD A 



Function 

Add to the A-address. The_A- 
address of the instruction (B) 
specified by B is increased by 
C. 



= 060 



ADD B 



Add to the B -address. The_B- 
address of the instruction (B) 
specified by B is increased by 
C. 



2 = 006 



ADD C 



Add to the C-address. The_C- 
address of the instruction (B) 
specified by B is increased by 
C. 



2 = 700 



SUB A 



°2 = 



070 



2 = 007 



SUB B 



SUB C 



Subtract from the A-address. 
The A-address of the instruction 
(B) specified by B is decreased 
by C. 

Subtract from the B -address. 
Analogous to SUB A, 

Subtract from the C-address. 
Analogous to SUB A. 



The sign of the instruction being modified remains unchanged and does not affect 
the outcome of the modification. Attempts to increase an address beyond 999 or 
decrease it below will result in erroneous operation not prevented by error stops. 

Three L -operations, TR A, TR B and TR C, complete the set of addres_s 
change operations. In each of them, the A-address specifies the instruction (A) 
to be called out and the B-address is the constant to which a specified address is 
to be compared. Incase of inequality, control is transferred to C. For example, 
if the instruction 6 750 325 200 (using L = 6, TR A) is given, control will be trans- 
ferred to 200 if the instruction in 750 reads 1 320 400 000 but control will proceed 
ahead if 750 contains 1 325 400 000. Summarizing: 



°1 = 



TR A 



Transfer on the A-address. The 
A-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 
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Numer. 



Alpha. 



Function 



Oj = 7 TR B Transfer on the B -address. The 

B-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 

°1 = 8 TR C Transfer on the C-address. The 

C-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 

As an introductory example, the summation in the section on LOOP OPERA- 
TIONS will be programmed again using address change methods. This would be 
an inefficient choice in an actual problem, but it will best illustrate the difference, 
as well as the analogy between the two methods. It is again assumed that register 
700 contains zero at the start, but the steps analogous to the resetting of the loop 
box will be included. 



Inst. 



Alpha. 



Sign 



0, 



A or 0. 



B 



101 


SET A 
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500 


103 


801 


102 


SET B 


+ 





050 


103 


901 


103 


MPY 


+ 


3 


[ ] 


[ ] 


000 


104 


ADD 


+ 

* 


1 


000 


700 


700 


105 


ADD A 


+ 





600 


103 


001 


106 


ADD B 


+ 





060 


103 


001 


107 


TR A 


+ 


6 


103 


821 


103 



108 



Next instruction in the problem. 



The brackets in the A- and B-addresses of instruction 103 are used to indicate that 
these addresses are variable and will be supplied by the program before the instruc- 
tion is executed, hence what is written there when the program is loaded into the 
machine is irrelevant. At the end of the program when instruction 108 is reached, 
memory location 103 will contain + 3 821 921 000. It is assumed that the summation 
just programmed is part of a larger problem in which it is used repeatedly. This 
is the reason for the SET A and SET B instructions. If 801 and 901 were simply 
loaded into their respective positions in instruction 103 initially, the summation 
would be performed correctly the first time it is used, but the next time when the 
summation is called for, instruction 103 would read + 3 821 921 000 and erroneous 
calculations would result. The SET instructions could, of course, have been in- 
serted after the completion of the summation, restoring instruction 103 to its proper 
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value for the next application, but this procedure is not recommended because it 
makes it more difficult to restart the problem from the beginning without reload- 
ing the program in case of an interruption (e. g. , error stop) during the loop. 

A more realistic example of the use of address change methods would be a 
calculation involving more than one summation index or parameter. Then, one 
of the fast and convenient LOOP orders would normally be used in the "inner loop", 
i. e. , the loop occurring most frequently, with address change operations con- 
trolling the "outer loop" or loops. Suppose, for example, that it is desired to 
calculate Sj ~ ™ a-x. for j = 1, 2, ... 5, where the a... are stored in 800 + 10^+ i 

i = l 
(i. e. , a n is in 811, a 12 in 812, etc. ; a in 821, a in 822, and so forth), the 
xj in 900 + i, and the S. are to be storedin 700+ j. It will be assumed that register 
500 contains zero. For completeness, the setting of all variable addresses to their 
initial values for repeated use of the summation program will be included. 



Instr. 


Alpha. 


Sign 


101 


SET A 


+ 


102 


SET C 


+ 


103 


MOVE 


+ 


104 


MPY 


_* 


105 


ADD 


+ 


106 


LOOP AB 


+ 



107 



110 



MOVE 



TR C 



+ 





9 



3 

1 




A or 0, 
c 

500 
005 
000 

I 1 

000 

110 
000 



107 



400 



706 



B 


C 


Comments 


104 


81 f 


Set variable 
' addresses to their 


107 


701 

J 


initial values 


500 


400 


Set register 400 
to zero 


901 


ooo"* 1 


"Inner loop" 


400 


400 


f i. e. , summation 


010 


104 


1 on i 



C 1 



108 


ADD A 


+ 





600 


104 


010 


109 


ADD C 


+ 





006 


107 


001 



103*1 



Store the result 

Increase address* 
for next repetitioi 
in the outer loop 
(j-loop) 

Test for end of 
j-loop 



111 Next instruction in the problem. 

A superficial examination of this program might suggest that only 1/5 of the 
program is devoted to actual arithmetic calculation (!), but it should be observed 
that in terms of the number of instructions executed by the machine when one 
complete summation is performed, the arithmetic ones are still in the majority, 
and in terms of execution time they comprise about 3/4 of the program. 
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In programming problems involving several loops, it may be helpful to 
consider the structure of a loop in terms of four phases: 

1. Initialization. Where addresses in the loop are set to their ini- 
tial values, registers used for summation are set to zero, etc. . 
The automatic resetting of the loop box and the fact that ad- 
dresses remain unchanged in memory tend to reduce the ini- 
tialization when the loop is controlled by a LOOP operation. 

In the summation problem above, steps 101 and 102 constitute 
the initialization for the outer loop, step 103 is the initialization 
for the inner loop. Notice that step 103 is repeated as a part of 
the outer loop. 

2. Execution. Comprising the mathematical operations of the loop, 
as well as any logical operations associated with a loop inside 
the one being executed. Above, the execution of the inner loop 
consists of steps 104 and 105 and the execution of the outer 
loop consists of 103-107. 

3. Modification. Where addresses, parameter values, etc. , are in- 
creased or decreased. The modification of the inner loop above 
is included in the LOOP instruction. The modification of the 
outer loop consists of steps 108 and 109. The position of the 
modification in the program in relation to the execution and test 
is frequently subject to choice. 

4. Test. Determining whether the loop is completed or further rep- 
etition is required. The LOOP instruction includes the test for 
the inner loop and step 110 is the test for the outer loop. 

Note: It is important to write loops in such a way that all ini- 
tializations are performed by the program, not by loading. If 
this rule is not followed, it will not be feasible to restart the pro- 
gram during testing or after a machine stop without reloading. 
For example, if a register is used for summation, it should be 
reset before being used in the summation loop by moving zero into 
it from another location, not by loading zero into it from a card. 

Many programmers find it helpful in programming a large problem to draw 
a block diagram or flow chart with one box representing each phase of each loop 
and arrows connecting the boxes showing the flow of control. 

The address change operations, particularly the SET operations, are fre- 
quently useful in non- repetitive situations as well. An example of this will be 
found in the section on SUBROUTINES. 

If a program appears to require a large amount of address modification 
and particularly, if this occurs because a quantity whose address is subject to 
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change is needed in many places in the execution of a loop, it may be advantageous 
to write the execution largely in terms of fixed addresses and perform the mod- 
ifications by moving data. Instruction 107, in the example above, illustrates this 
in a simple way: If the registers 700+ j themselves had been used in the summation 
process, (step 105) both the B- and C-addresses of step 105 would have required 
modification in the outer loop, as well as the C-address of instruction 103. For 
cases where several numbers are to be moved at the same time, a more general 
MOVE operation than the MOVE 000 used so far is available and will be described 
in the next section. 

IV. 4. MOVE 

The MOVE operation is defined as follows: 



Numer. 



0l = 9 



Alpha. 
MOVE 



Function 

Move. If A 4 s 000, the block of 
A consecutive words beginning at 
B is moved into the set of A con- 
secutive locations beginning at 
C. The words in the original 
locations are not destroyed, 
except where the two regions 
overlap. The number in location 
000 ("previous result") is not 
affected when C £ 000. Both 
C > B and C < B are permissible. 
An error stop occurs if C 4- A — 
12:1000. If A= 000, the word 
(B) specified by the B -address is 
moved into location C and into 000. 
It also remains in location B. 



MOVE With A = 000 differs from MOVE with A = 001 only in that the execution time 
with A= 000 is shorter and that the previous result location is affected. 

Note: If a number is to be moved from location B into 000 for use in a TR 
SGN or TR EXP operation on the next step, MOVE 9 001 B 000 must not be used, 
since these transfer operations work strictly according to their definitions (see 
Sec IV- 1. )• The correct instruction would be MOVE 9 000 B 000. (Internally, 
these transfer operations inspect a duplicate "previous result" location rather than 
000!) 

As an example, suppose Xj is in 701, x 2 in 702, . . . , x 5 in 705 and the instruc. 
tion MOVE 9 005 701 703 is given. Then x x will be found in 703, x 2 in 704, .... 
x 5 in 7 07, after execution. 

In conclusion, it should be pointed out that the use of the logical operations 
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is by no means restricted to the straightforward functions for which they are 
primarily intended. The programmer will find innumerable ways of increasing 
the efficiency and elegance of his programs by unusual applications, particularly 
of the address change operations. As a weird example, suppose it is desired 
to multiply the numbers located in registers 1, 4, 9, 16, 25, 36, ..., 400(f) 
by a constant located in 600 and store the results in 501, 502, 503, . .'. , 520: 

Instr. Alpha. Sign 0j A or 2 B C 

898 SETA + 500 9 00 001 

899 SETC + o 005 901 003 
I™ MPY - 3 [ ] 600 501 

901 ADDA + 600 9 00 T "I 

902 ADD C + 006 9 01 002 

903 LOOP C + 001 020 900 

904 Next instruction in the problem. 
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V INPUT-OUTPUT OPERATIONS 



V. 1. CARD FORM 

By a card form is meant a specific assignment of card columns to form fields 
for data, instructions, identification, etc. , in connection with a given program or 
interpretive system. In the 650, information is transmitted to and from cards 
through a control panel, and anyone whose needs call for a special card form can 
adapt it for use in connection with the interpretive system merely by simple control- 
panel wiring. For most needs, the following card form, associated with the inter- 
pretive system, is likely to be found adequate. At this point, only brief definitions 
of the card fields will be given for reference in subsequent sections where their 
use will be explained in detail: 



Columns 



Definition 



1-4 
5-6 

7-9 

10 

11 
12-21 

22 
23-32 

33 
34-43 

44 
45-54 

55 
56-65 

66 
67-76 
77-79 

80 



Card number 

Deck number 

Location 

Word count 

Sign of word 1 

Word 1 

Sign of word 2 

Word 2 

Sign of word 3 

Word 3 

Sign of word 4 

Word 4 

Sign of word 5 

Word 5 

Sign of word 6 

Word 6 

Problem number 

Tracing identification 



The same card form is used in all input-output operations, as well as in trac- 
ing. Both instructions and data are signed ten-digit words and are entirely indistin- 
guishable in connection with input-output operations. 

V. 2. PUNCHING 

At any point in the problem, the machine may be ordered to punch into cards 
the contents of any set of memory locations, together with appropriate identification. 
In some problems, it may be desirable to punch out answers one at a time, perhaps 
together with the values of relevant parameters; in others it may be preferable to 
punch out a large amount of information at less frequent intervals. There are also 
cases where it is advantageous to punch out instructions: In connection with testing 
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(see PROGRAM TESTING) in order to examine a program interrupted at a chosen 
point, and in connection with loading, (see LOADING) in order to reduce the size of 
a deck of cards. All of these ends are served by the following instruction: 

°2 = 41 ° PCH Punch cards. The block of consecutive 

words beginning at B arid ending at C 
(inclusive) is punched into cards. Five 
words and a word count of 5 are punched 
into each card but the last, whose word 
count will be the remainder when C-B + 
1 is divided by 5. On each card, the 
location from which word 1 was punched 
is punched into columns 7-9. The 
words in storage are not destroyed. A 
cumulative count of the number of cards 
punched during the problem (i. e. , since 
LOADING) is punched into columns 1-4. 
The problem number (see LOADING) is 
punched into columns 77-79 and zero is 
punched into columns 6 and 80. An error 
stop occurs if B > C. 

If it is desired to punch six words to a card, this may be done by adding a 
special card behind the punching deck (see LOADING). This card should have an 
x-punch in column 5, 1969 in columns 6-9, 1 in column 10, a 12 punch in column 
11, and 00 0006 0000 in column^ 2-21. 

The punched cards are likely to be used for one (or both) of two purposes (in 
addition to possible processing on other equipment): The information on them may 
be printed on a tabulator or they may be loaded (or READ) into the 650 at a later 
time. Details of the printing will not be given here, since they depend on the char- 
acteristics of the tabulator, but the printing form may be assumed to be roughly 
identical with the card form with proper spacing between words. (Suggestions on 
tabulator wiring are given in Sec. XI. 5. ) It is assumed that the suppression of the 
superfluous words punched into the last card, if its word count is not 5, will be 
performed on the tabulator control panel. If this is not feasible, it may be done in 
the 650 by adding three cards to the punching deck. For details, see Sec. XI. 1. 

Selective spacing between lines in printing may be accomplished in several 
ways, even though no operation in the 650 is provided for this purpose. A brief 
discussion will be given here, since spacing considerations may affect the use of 
the PCH operation in programming. Through the setting of switches on the tab- 
ulator, a choice of any of the following spacing alternatives may be provided: 

(a) Single or double spacing. 
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(b) Spacing between every n lines (with n chosen by wiring, nor- 
mally, e.g., n = 10). 

(c) Spacing after any line whose word count is less than the word 
count of the preceding line. 

(d) Spacing before any line whose location number has a units digit 
smaller than the units digit of the location number of the pre- 
ceding line. 

Alternative (c) is suited for the printing of information punched from fairly 
large blocks of locations by one PCH order. Spacing will occur after each block, 
unless the block length is a multiple of 5, which can be avoided by programming. 
Alternative (d) is intended for information punched repeatedly from the same set 
of locations and provides the option of spacing when the loop is interrupted, e. g 
for changing a parameter value. 

V. 3. LOADING 

When a program has been written, and careful inspection reveals no further 
errors, it is key punched into cards following the card form given in Sec. V. 1. To 
reduce to a minimum the number of errors to be found with the aid of the 650, the 
cards should be run through a verifier operated by another person or, alternatively, 
key punched independently by two operators and compared on a reproducer. The 
programmer has the option of specifying the number of words to be punched to a 
card: Punching 5 or 6 to a card will keep the program deck small from the outset 
and eliminate the need for condensing the deck on the 650 later. Punching one word 
to a card is felt by some programmers to facilitate changes. Each card must have 
in columns 7-9 the location into which word 1 is to be loaded, and in column 10 the 
number of words to be loaded from the card into consecutive locations. Columns 
1-6 and 77-80 are not read by the 650 (except that the problem number is read from 
the last card, see below) and may be used by the programmer as he deems best. 
Each column of each field to be used by the machine must contain one and only one 
punch and an error stop is provided to enforce this rule. A 12-punch is used for 
plus, an 11- or x-punch for minus and a 0-punch — not a blank column- -for zero. If 
the word count is less than 6, unused word fields and sign columns may be left 
blank. No distinction is made between data and instructions in key punching and load- 
ing. 

LOADING is the process of feeding data and instructions into the machine at 
the beginning of a problem. If the previous user of the 650 was not using the inter- 
pretive system or if there is any reason to doubt that the system is correctly 
stored on the drum, the program deck should be preceded in loading by a deck which 
loads the interpretive system (in 51. 9 seconds) into the memory locations aboVe 999. 
Before the program deck, the programmer may also place a Reset Memory Card, 
which will (in 6. 3 seconds) reset each of the memory locations 001-999 to minus 
zero. (This is useful in connection with the punching out of sections of memory in 
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testing. ) Immediately behind the program deck--no blank cards are used in the 
card reader in connection with this interpretive system- -the programmer places 
one of two nine-card decks to inform the machine whether he wants normal opera- 
tion or TRACING described in a later section. (If he knows that he wants the 
same mode of operation as the previous user, he can omit these cards but the 
gain is only 2. 7 seconds. ) Last, he must place a so-called transfer card with a 
zero punched in column 10, the problem number in 77-79 and the location of the 
instruction at which the program begins in columns 7-9. The word fields on this 
card may be left blank. 

The loading program automatically resets the loop box, the card counter 
(see PUNCHING) and location 000 to zero. 

The order in which the program cards are loaded is irrelevant, unless the 
same location is loaded into from more than one card, in which case the last such 
card, of course, determines the contents of the location. This may occur in 
connection with changes of a temporary nature, which may be placed at the end of 
the deck and later removed, leaving the program in its original form. In the deck 
which loads the interpretive system, the order of the cards must be preserved, 
and an error stop is provided to insure this, thereby ascertaining that no part of 
the system is missing. 

In summary, a complete deck to be loaded must contain: 

System deck (173 cards) 
Reset Memory card (optional) 
Program deck 

Mode-of-operation deck (9 cards) 
Transfer card 

The control console of the 650 need be of almost no concern to the user of 
the interpretive system under normal conditions. He must only make sure that 
all switches on the console are set in a fixed manner required by the system, and 
these settings will now be listed without any description of the function of the 
switches. Certain ways of using the console are described in the sections on READ- 
ING and PROGRAM TESTING. 

Switches Settings 

Storage entry 70 1951 1333 -f 

Programmed stop Stop (see COND STOP) 

Half cycle Run 

Address selection 1338 (see STOPS) 

Control Run 

Display Upper Accumulator 

Overflow Stop 

Error Stop 

To start a problem, the deck to be loaded is placed in the card reader, and 
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the following keys are depressed in order: 

(1) COMPUTER RESET 
(2) PROGRAM START 
(3) START 



(on the console) 
(on the console) 
(on the card reader) 



When the last card leaves the hopper, the machine stops and the key labelled 



(4) END OF FILE 



(on the card reader) 



is depressed. If the deck has been correctly put together, the execution of the pro- 
gram will then start automatically. 

The program deck may be run out at any time after loading by depressing the 
START key, unless a READ instruction is contained in the program. Blank cards 
should be inserted into the PUNCH hopper and the START key on the punch side 
depressed. 

To make the 650 produce a condensed program deck in case the program was 
originally key punched one instruction to a card, a PCH instruction should be given 
at the very beginning of the program. This instruction may be bypassed during 
subsequent executions of the program merely by changing the location number on 
the transfer card. 

V.4. READING 



In some problems, particularly in applications of a data processing nature, 
it may be desirable to read information into the machine during the execution of 
the program without manual interference. This is accomplished by the READ 
operation: 



2 = 400 



READ 



Read cards. The block of consec- 
utive storage locations beginning 
at B and ending at C (inclusive) 
is read into from cards. The ad- 
dress B must appear in the location 
field on the first card, as well as 
in the READ instruction, and the 
location field on each card follow- 
ing must contain the sum of the 
word count and location on the 
previous card. The sum of the 
word counts of all cards to be 
read must be C-B+ 1. Violations 
of these requirements, which have 
been included for the programmer's 
protection, will result in error 
stops. 



-29- 



The cards to be read should be placed in the hopper of the card reader imme- 
diately following the transfer card (no blank cards). 

The decisions made with the aid of conditional transfers and other logical 
operations are normally based on criteria predetermined by the programmer and 
incorporated into the program. If the programmer wishes to influence the pro- 
gram during its execution, e. g. , on the basis of a result displayed on the console 
in connection with a COND STOP instruction, he may do so using the CONS opera- 
tion: 

C>2 — 401 CONS Read console. The machine stops 

if the PROGRAMMED STOP switch 
is on STOP. Zero is displayed on 
the display lights and 1131 on the 
address lights. When the PROGRAM 
START key is depressed, the num- 
ber entered on the STORAGE 
ENTRY SWITCHES is stored in 
location C and in 000 (the "previous 
result" location). If the PRO- 
GRAMMED STOP switch is on 
RUN, the storing takes place with- 
out a stop preceding. The B -ad- 
dress is ignored but should be 
filled, e. g. , with zeros. 

An example of an application of CONS might be the feeding in of an "educated 
guess" for a starting value in connection with the solution of algebraic equations. 
Another application, involving only the storage entry SIGN switch, might be to 
continue a program until another user is ready to take the machine, at which time 
a change in the SIGN switch setting, interpreted by a TRSGN operation, causes the 
program to punch out intermediate results for later restart. 
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VI PROGRAM TESTING 



VI. 1. MEMORY PRINT -OUT 

The choice of methods for testing ("debugging") a program by comparing 
results of machine calculation to known quantities or to results of independent 
calculations by other means is governed by the relative availability of machine 
time and programmer time. If machine time is freely available, testing with the 
aid of the control console is highly efficient, as well as instructive and enjoyable, 
as soon as a certain facility for operating the console has been acquired. Par- 
ticularly in the case of small problems, the method of tracing- -where a card is 
punched for each instruction executed, showing all numerical and logical quantities 
associated with the execution- -may be the most desirable in that it gives an almost 
certain clue to the difficulty within a predictable, if not very short, period of 
machine time and allows the programmer to study the material at his leisure. 
The method most economical of machine time and yet frequently sufficiently 
illuminating is that of memory print-out. It might be suggested that on most prob- 
blems in a busy but not heavily over-leaded installation, the methods be used in 
the order reverse to that in which they were mentioned here. Some directions 
for their use will now be given. 

The memory print-out method simply consists of inserting temporarily into 
the program at one or several suitably chosen points PCH orders (see PUNCHING) 
calling for the punching of blocks of information- -data or instructions --which, when 
printed on the tabulator, will give a picture of the progress of the program. Since 
1000 words may be punched 5 to a card in two minutes, it is not out of the question 
to punch out the contents of every register used in a problem- -including all the 
instructions — several times. To get the most benefit from this method, the pro- 
grammer should, in any problem that does not threaten to fill the entire available 
memory, avoid using the same storage location for storing different quantities at 
different times whenever feasible, so that as many partial results as possible are 
preserved for the memory print-outs. Whenever a test case of a problem is run, 
even if memory print-out is not chosen as the primary testing method, it would 
certainly be advisable to make the last instruction of the test deck punch out the 
entire memory used. A flexible alternative would be to have scattered through the 
program CONS--TR SGN combinations which transfer control to a PCH order if 
the storage entry sign switch is turned to minus. 

Temporary instructions may be inserted into a program in two ways: Either 
they are included in the normal sequence of instructions when the program is 
initially written and replaced either by NOOP instructions (see Sec. IV. 1) or by 
transfer to the next non-temporary instruction when no longer needed,or else one 
of the regular instructions of the program is replaced by a TR to a vacant location 
L, the regular instruction is placed in L, the temporary ones inL+1, L+ 2, etc. , 
and at the end of this temporary sequence a TR back to the normal program is given. 
In either case, the temporary instructions may (as suggested in LOADING) be kept 
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Colum: 


ns 


1-4 




6 




7-9 




10 




11-21 





as a separate deck at the end of the program deck, eliminating any changes in the 
main program deck and simplifying bookkeeping. 

VI. 2. TRACING 

If the tracing deck of nine cards is loaded with the program deck, (see LOAD- 
ING) the machine will automatically start tracing from the beginning of the program, 
as specified by the transfer card. Before the execution of each instruction, a card 
with the following information will be punched: 

Definition 

Card number (cumulative) 

Zero 

Location of the instruction 

about to be executed. 

Six 

The instruction as stored in 

memory. 

22-32 The instruction as modified 

for execution (i.e., with the 
loop box added if minus). 

33-43 The contents of the loop box. 

44-54 A if A * 000, zero if A =000. 

55-65 B 

66-76 The contents of location 000 

(i. e. , the result of the last 
mathematical MOVE 000 or 
CONS operation). 

77-79 Problem number. 

80 Eight (used by the tabulator 

for automatic selection of a 
different printing form for 
trace cards). 

The punching rate will be 100 cards per minute except in the case of very time-con- 
suming operations, such as, the moving of a large block of information. The ad- 
vantage of punching the trace card before execution is that information will be 
punched for an instruction whose execution is interrupted by an error stop. In the 
case of instructions (such as LOOP or TR EXP) whose B-address does not refer 
to a memory location, the quantity B is irrelevant. Tabulator wiring to suppress 
the printing of B in such cases can be provided if sufficient selection equipment is 
available. The PCH operation is bypassed when the machine is operating in the 
tracing mode, i. e. , PCH is equivalent to NOOP. 

If a program is too long to be traced in its entirety or if this is unnecessary, 
selective tracing may be effected by using the following operations: 
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2 = 450 



START TR 



o = 451 
2 



STOP TR 



°2 = 452 



ST TR ERAS 



Start tracing. If the nine -card 
tracing deck has been loaded, 
the machine will start tracing 
from the next instruction. ..If it 
is already tracing, it will con- 
tinue to trace. The B and C 
addresses are ignored. If the 
deck for normal operation has 
been loaded, START TR will 
be equivalent to NOOP. 

Stop tracing. If the machine is 
tracing, it will discontinue trac- 
ing immediately. If it is not 
tracing, STOP TR will be equiv- 
alent to NOOP. The B and C 
addresses are ignored. 

Start tracing and erase itself. 
If the tracing deck has been 
loaded, the machine will start 
tracing from the next instruc- 
tion. If it is already tracing, 
it will continue to trace. If the 
deck for normal operation has 
been loaded, . tracing will not 
begin. In all cases, the ST TR 
ERAS instruction will be replaced 
in memory by a NOOP (0 2 = 454) 
during its first (and only!) execu- 
tion. The B and C addresses are 
ignored. 



The bypassing of the PCH operation is in effect as long as the trace program is on 
the drum and is not affected by the selective tracing orders. To make PCH opera- 
tive, the nine-card deck for normal operation must be loaded. 

The purpose of the ST TR ERAS operation is to make it possible to trace the 
repetitive steps of a loop either once or twice and then stop tracing until the loop 
is completed. To get the steps traced once, one may place the pair STOP TR, ST 
TR ERAS at the beginning of the repeated portion of the loop; to get them traced 
twice, one places this pair of instructions at the end immediately preceding the 
test. As a specific example, suppose it is required to trace twice the steps of the 
loop programmed in the section on LOOP OPERATIONS and suppose vacant loca- 
tions are available from 900 up. Assume that the machine is tracing as it enters 
the loop. The original program reads as follows: 
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101 
102 
103 

104 



MPY — 3 

ADD "+" 1 

LOOP AB + 

Next instruction in the problem. 



801 


901 


000 


000 


700 


700 


110 


020 


101 



The following instructions could be added as a temporary deck at the end of the pro- 
gram deck: 



102 


TR 


+ 





203 


000 


900 


900 


ADD 


+ 


1 


000 


700 


700 


901 


STOP TR 


+ 





451 


000 


000 


902 


ST TR ERAS 


+ 





452 


000 


000 


903 


TR 


+ 





203 


000 


103 



Notice that the TR instruction gets loaded into 102 after the regular program, replac- 
ing the ADD instruction, as explained in LOADING. This example is, of course, 
unrealistic in that selective tracing would hardly be needed for testing such a simple 
loop. 

VI. 3. CONSOLE TESTING 

Testing with the aid of the control console requires some familiarity with the 
internal structure of the interpretive system (see Sec. XI) and with the basic language 
of the 650. Console testing is more attractive on the 650 than on most machines due 
to the ADDRESS STOP feature: If the CONTROL switch is turned to the ADDRESS 
STOP position, the execution of the program will proceed at electronic speed until 
the address set up on the ADDRESS SELECTION switches is reached. At that point, 
the machine stops, and the contents of various registers may be examined on the 
display lights or the program may be continued manually one step at a time. Alter- 
natively, the program may be punched out on cards at this point by merely feeding 
in one card with a PCH instruction, going into any vacant location, followed by a 
transfer card specifying this location. Console testing, in connection with the inter- 
pretive system, is likely to be needed only in exceptional cases. 

The ADDRESS STOP feature of the 650 may be used in conjunction with a spe- 
cial address stop transfer card when it is desired to start tracing from a certain 
instruction N in the middle of a program after running at full speed up to that point. 
(This may, of course, alternatively be accomplished using the tracing operations 
described in Sec. VI. 2, but then the value of N must be decided upon in advance 
and the proper program changes key punched. ) The procedure is as follows: Set 
the ADDRESS SELECTION switches to N and turn the CONTROL switch to ADDRESS 
STOP. Load as usual and run until the machine stops at the instruction N. (For 
details on possible earlier stops see below. ) Then set the CONTROL switch to 
RUN and load the tracing deck followed by the address stop transfer card. Tracing 
will begin immediately and the first instruction traced will be N. 

In choosing N it must be remembered that the loop box and location 000 are 
reset to zero when the tracing deck is loaded. If this restriction is inconvenient, 
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it can be circumvented by placing a special card in front of the tracing deck. The 
card counter and the problem number are also reset to zero, unless the tracing 
deck has been modified to prevent it. 

If the CONTROL switch is kept in the ADDRESS STOP position when the pro- 
gram deck is loaded, one stop will occur when location N is reset by the memory 
reset card and another when the programmer's instruction is loaded into N. Also, 
stops may occur before instruction N is reached in the program, if N is referred 
to in an ADDRESS CHANGE or MOVE operation (but not if N is one address in a 
conditional transfer instruction and control is transferred to the other address). 
After each stop, operation will resume when the PROGRAM START key is depressed. 
If the CONTROL switch is left in the ADDRESS STOP position during tracing, two 
stops will occur each time N is referred to (and one if N is the B-address of a 
transfer instruction). 

The program can be continued at full speed (punching mode) after a period of 
tracing by following the procedure described above with the punching deck in place 
of the tracing deck. 

The address stop transfer card has 69 1976 1952 24 1061 1098 in columns 
1-20 and a 12-punch in each of columns 1, 10 and 20. The special card for by- 
passing the resetting steps in loading has. 69 1953 1952 24 1278 1953 70 1951 1344 
in columns 1-30 and a 12-punch in each of columns 1, 10, 20 and 30. (See Deck 7, 
Sec. XI. 7.) 

If the value of N has been decided upon in time to get it key punched into a 
regular transfer card, (Sec. V. 3) this card may, of course, be used in place of 
the address stop transfer card in the procedure described above. 

A programmer familiar with the internal structure of the interpretive system 
will be able to think of many other cases where special needs can be met using 
maching^ language cards ("load cards"). 
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VII SUMMARY OF OPERATION CODES 



1 OPERATIONS 



2 OPERATIONS 



Num. 



1 
2 
3 
4 
5 



6 
7 
8 



Alpha. 

GO to 

ADD 

SUB 

MPY 

DIV 

NGMPY 

TR A 
TR B 
TR C 

MOVE 



Page 
Ref. 

8 

10 
10 
10 
10 
10 

19 
20 
20 

23 



Num. 

000 
200 
201 
202 
203 
204 
205 

100 
010 
001 
110 
101 
011 
111 

500 
050 
005 
600 
060 
006 
700 
070 
007 





Page 






Page 


Alpha. 


Ref. 


Num. 


Alpha. 


Ref. 


UNC STOP 


13 


300 


SQRT 




COND STOP 13 


301 


EXP E 




TR SGN 


14 


302 


LOG E 




TR EXP 


14 


303 


SIN R 




TR 


14 


304 


COS R 




TR SUBR 


15 


305 


ART R 




TR OUT 


15 


350 
351 


ABS 
EXP 10 




LOOP A 


16 


352 


LOG 10 




LOOP B 


16 


353 


SIN D 




LOOP C 


16 


354 


COS D 




LOOP AB 


17 


355 


ART D 




LOOP AC 


17 








LOOP BC 


17 


400 


READ 


29 


LOOP ABC 


17 


401 


CONS 


30 






410 


PCH 


26 


SET A 


18 








SET B 


18 


450 


START TR 


33 


SET C 


18 


451 


STOP TR 


33 


ADD A 


19 


452 


ST TR ERAS 33 


ADD B 


19 








ADD C 


19 


454 


NOOP 


15 


SUB A 


19 








SUB B 


19 








SUB C 


19 
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VIII STOPS 

Error circuits in the 650 will stop the machine on attempts to use invalid 
information, such as, that represented by blank columns or double punches, as 
well as on several kinds of machine malfunctioning, and will indicate on the control 
console the nature of the error. If this occurs during the loading of a new deck, 
the cards should be examined. In other cases, a note should be made of the indica- 
tions on the console, and the procedure that led to the stop should, if possible, be 
repeated exactly in order to determine whether the error is systematic m nature. 

All stops, which are part of the interpretive system, will now be listed. 
Conditional stops will occur only if the PROGRAMMED STOP switch is set to 
STOP. On a conditional stop, the PROGRAM LIGHT in the OPERATING section 
of the console will be on and no lights in the CHECKING section should be on. The 
program will continue if the PROGRAM START key is depressed. On an uncon- 
ditional stop, the STORAGE SELECTION light in the CHECKING section will be on. 
Normally, operation should be discontinued after an unconditional stop and changes 
made in the program in order to avoid the stop. Alternatively, the program may 
be continued by having a transfer card (see LOADING) in the card reader, specifying 
the instruction to which control should proceed when the COMPUTER RESET and 
PROGRAM START keys are depressed. 

The location of the interpretive system instruction xxx on which the machine 
has stopped, may be determined by displaying the contents of location 1098 on the 
console. The display lights will show 60 Oxxx 1107. This process, called "monitor- 
ing", may be performed as described in the 650 manual or, alternatively, by setting 
the storage entry switches to 60 Oxxx 8000 and depressing the COMPUTER RESET, 
PROGRAM START and PROGRAM STOP keys. 

If, in an exceptional case, it would be advisable to proceed to the next instruc- 
tion after an unexpected unconditional stop, this may be done manually as follows: 

(1) Set the CONTROL switch to MANUAL. 

(2) Check that the ADDRESS SELECTION switches are set to 1338. 

(3) Depress the COMPUTER RESET key. 

(4) Depress the TRANSFER key. 

(5) Set the CONTROL switch to RUN. 

(6) Depress the PROGRAM START key. 

As a result of this procedure, zero will be stored at C and 000 before the next instruc- 
tion is executed. If this is not desired, the ADDRESS SELECTION switches should 
be set to 1095 in step (2). To repeat the same instruction (on which the stop occurred) 
the switches are set to 1098. 

There is an alternative manual procedure for restarting after an unconditional 
stop which is simpler in the case of frequent use but is not recommended in general 
because it requires changing the setting of the STORAGE ENTRY switches. They 
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are used in LOADING and must be set back to their normal positions for the next 



user 



(1) Set the STORAGE ENTRY SWITCHES to 00 1951 1338-1- (or 
00 1951 1095+ if zero is not to be stored or 00 1951 1098 + 
to repeat). 

(2) Depress the COMPUTER RESET key. 

(3) Depress the PROGRAM START key. 
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CONDITIONAL STOPS 
Address Lights 

1120 

1131 

1715 

1835 



Normal Cause 

Programmed COND STOP. (Display lights show B. ) 

CONS (Check STORAGE ENTRY switch setting. ) 

Loss of accuracy in SIN (Exponent of B exceeds 52) 
or COS. 

Loss of two or more digits of accuracy in LOG. 



UNCONDITIONAL STOPS 



2222 



(MOVE with C + A - 1 >_ 1 000. 
JpCH with B > C + 1. 

READ with incorrect loc. or word count. 



3333 
4444 

5555 



6666 



7777 



9999 



DIV with B = 0. 
SQRT with B < 0. 

MPY with result out of range. 
DIV with result out of range. 

^EXP with result out of range. 
LOG with B <_0. 

SIN with exp. of B exceeding 58. 
COS with exp. of B exceeding 58. 

Cards missing or out of order in the system deck 
being loaded. 

Programmed UNC STOP (Display lights show B). 
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IX EXECUTION TIMES 



a ?% eX ,fn C r i0n timCS HSted in tMs SeCtion are based on the standard 650 drum 
speed of 12, 500 r. p. m. They represent approximate theoretical estimates derived 
in the case of the mathematical operations, from simple assumptions regarding the 
distribution of the numbers to be operated upon. For example, the part A, of a 
floating -decimal number A = A, . 10 1 is assumed to be uniformly distributed 
between 1 and 10, although in physical problems there are reasons that favor a 
logarithmic distribution; extremely small and extremely large exponents are con- 
sidered very unlikely, etc. It is further assumed that the programmer has chosen 
storage locations on the drum without regard to timing, ignoring the fact that in 
the case of some operations the execution time will be a few milliseconds shorter 
for numbers stored in certain sections of memory. Some, but not nearly all, of 
tne time estimates have been verified by tests. 

It should be stressed that the estimates of execution times are needed only 
for making comparisons or estimates of running time for problems or for choosing 
efficient ways of programming and will never affect the result of an operation. In 
comparing these estimates to estimates given for other interpretive systems or 
subroutines it is important to verify by sample calculations or machine tests that 
the assumptions are realistic. 



r« ^ minimize the Size of the table > the execution times listed refer to a basic 
table COrreCtl ° ns to be add *d in other cases are given at the beginning of the 
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650 INTERPRETIVE SYSTEM 
ESTIMATED AVERAGE EXECUTION TIMES IN MILLISECONDS 



(a) If A^ 000, add 7. 2 ms. for ADD and SUB, 6. 3 ms. for MPY, NGMPY and 
DIV*. 

(b) If C £ 000, add 6. 1 ms. for all mathematical operations, MOVE 000 and 
CONS*. 

(c) If the instruction has a minus sign, add 4. 8 ms. for all operations. 

(d) If, after a TR EXP or LOOP operation, control will proceed to the next 
instruction rather than to C, add 4. 8 ms. 



ADD 


65.7 


UNC STOP 


28. 8 


SQRT 


206 


SUB 


65.7 


COND STOP 


29.8 


EXP E 


197 


MPY 


67.2 


TR SGN 


19.2 


LOG E 


202 


DIV 


74.3 


TR EXP 


24. 


SIN R 


192 


NGMPY 


67.2 


TR 


19. 2 


COS R 


187 






TR SUBR 


44.4 


ART R 


238 






TR OUT 


26. 


ABS 
EXP 10 


33.2 
187 


TR A 


37. 3 


LOOP A 


24. 


LOG 10 


207 


TRB 


37. 3 


LOOP B 


28. 8 


SIN D 


240 


TR C 


42. 1 


LOOP C 


24.0 


COS D 


235 






LOOP AB 


24. 


ART D 


271 






LOOP AC 


24. 






MOVE 00037.7 


LOOP BC 


28.8 


READ 


One card: 


MOVE 


40. 8 + 12A 
(A = no. of 


LOOP ABC 


24. 




101-4- 14n 

(n= no. of words. ) 




words. ) 


SET A 


55. 3 




Succeeding cards: 






SET B 


55. 3 




300 each. 






SET C 


55. 3 


CONS 


28. 8 






ADD A 


44.9 


PCH 


One card: 






ADD B 


44.9 




163+ 12. 5n 






ADD C 


44.9 




(n = no. of words. ) 






SUB A 


44.9 




Succeeding cards: 






SUB B 


44.9 




600 each. 






SUB C 


44.9 


START TR 
STOP TR 
ST TR ERAS 
NOOP 

TRACING 

LOADING 


28. 8 
24. 
38.9 
24.0 

600 per card. 

300 per card. 



*(See next page for footnote. ) 
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*Thbse who are particularly interested in time considerations may wish to know 
the exact increments on which the weighted averages in (a) and (b) are based: 

In ADD and SUB, 4.8 ms. if 17 < A< 41 (mod 50) 

9.6 ms. if 1 < A< 16 or 42 < A< 49 (mod 50) 

In MPY, NGMPY and DIV, 4. 8 ms. if A > 17 or A = 1 (mod 50) 

9. 6 ms. if 2 < A< 16 (mod 50) 

In all mathematical operations, MOVE 000 and CONS: 

4. 8 ms. if 7 < C < 42 (mod 50) 

9. 6 ms. if 1 < C < 6 or 43 < C < 49 (mod 50) 

An easily remembered programming rule could be extracted from this in- 
formation: If locations between 17 and 41 (mod 50) are used for storing numbers, 
the increments given in (a) and (b) may be replaced by 4. 8 ms. 



li 
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X SPECIAL, TOPICS IN PROGRAMMING 



X. 1. SUBROUTINES; TRANSLATION 

A subroutine is a program expected to be of use as a part of the program in 
several problems or in several sections of the same problem. The mathematical 
operations in the interpretive system are indeed subroutines written in the basic 
language of the machine and reached through their operation codes, and anyone 
desirous of preparing an additional subroutine of this type may avail himself of a 
vacant operation code (see Sec. XI) and write the program in a part of the memory 
below 1000. 

Subroutines written wholly or partly in the interpretive language may be 
reached conveniently using the TR SUBR operation defined in Sec. IV. 1. Suppose 
the subroutine begins at 900 and ends at 9 35. Instruction 900 should read: SET 
C 005 935 C 3and instruction 935 should read: TR 203 000C 1 Now suppose 
the subroutine is needed at step 700 in a program, and when it has been used, 
control is to be transferred to 680. Instruction 700 should read: TR SUBR 204 
680 900. The TR SUBR operation will take the quantity ("return address") 680, 
place it in the C-address of instruction 900 and then transfer control to 900. Instruc- 
tion 900, in turn, places 680 in the C-address of 935, and when instruction 935 is 
reached at the end of the subroutine, it transfers control to 680 as originally desired. 
Hence, the programmer using the subroutine only needs to know the identifying 
number 900; the transfer of control to and from the subroutine is handled by the 
TR SUBR in conjunction with the two instructions 900 and 935 provided by the writer 
of the subroutine. Subroutines needing only one input number and giving only one 
result (such as, the evaluation of one Bessel function for a given value of the argu- 
ment) will normally assume the input to be in 000 and will deliver the result there; 
in the case of several numbers, specified locations normally within thejblock 
occupied by the subroutine would be used forTnput"and/or results. Subroutines may, 
of course, be used inside other subroutines without restriction. 

If the locations occupied by a subroutine are needed for another purpose, e. g. , 
another subroutine in the same problem, the subroutine may be translated to a 
different set of locations by a translating program developed by Miss D. C. Leagus. 
When the subroutine is written entirely in the interpretive system, the programmer 
is required only to separate data and constants from instructions, and the translating 
program will automatically decide which addresses of each instruction are subject 
to translation. Machine language instructions may also be used in a subroutine to 
be translated, provided certain conditions specified by the translating program are 
adhered to. 

Subroutines for the solution of cubic equations and of systems of linear 
equations have been written at the Laboratories. 
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X. 2. UNNORMALIZED INPUT; TRANSITION BETWEEN FLOATING- AND FIXED- 
DECIMAL FORM 

Nearly all of the mathematical operations in the system assume that the float- 
ing decimal numbers to be operated upon are in the normalized form defined in 
Sec. II. 2, i. e. , that the leading digit is different from zero unless the entire num- 
ber is zero. In processing empirical data, key punching is often facilitated by 
permitting leading zeros and reproducing a constant exponent. Such unnormalized 
data may be used in the interpretive system provided the first operation in which 
it is used is ADD or SUB with operand exponents differing by less than 10. 

A special case of unnormalized input is that of a zero with a non-zero machine 
exponent. If such a zero is added to a non-zero number with a smaller exponent, 
a number of digits equal to the difference between the exponents are lost. Con- 
sequently, zero should be equipped with exponent 00 unless the programmer knows 
in detail how the zero will be used in his program. Special provisions in the MPY 
and DIV routines make it possible to use a zero with machine exponent 00 in tnem 
without danger of exceeding the exponent range negatively. 

The converse problem of producing unnormalized output, e. g. , for the print- 
ing of tables in fixed-decimal form or for calculations in machine language is easily 
solved at the expense of one digit. Suppose for example that the numbers N. to be 
"unnormalized" or "fixed" are known to be less than 10 4 and output in the form 
XXXX.xxx is desired. Add the number 1000000054 (i. e. , 10, 000. ooo) to N. if 
Ni > 0, subtract it from N. if N. < (using TR SGN) and punch. The output of the 
f orm + (10, 000+ |Nj|} , is ready to be printed on the tabulator with the leading 
1 and the constant exponent 54 suppressed by hammerlocks or wiring. If the num- 
bers are to be used in machine language, the 1 and 54 are shifted out. Rounding 
to a smaller number of digits is obtained by choosing the exponent of the additive 
constant (1000000054) correspondingly larger. 

X. 3. EXAMPLES 

In conclusion, two problems will be programmed in order to illustrate the use 
of many of the operations and methods described. 

First, suppose it is desired to evaluate the "error function", 

(i) $Cx) = -^. f e- t2 <Lt 

for a set of values x = a, a+& , a+2A, . . . , a+10A, using the RAND approximation 

(2) <£*(x>= i -Ca 1 n+a E n 2 +a 3 n 3 -+-^n 4 +a 5 ri ? )^ , cx) ) 
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where 



(3) 



ft- 1/O+px)^ (pis 3 numerical constant). 



and 



(4) Q.'CX)=^i = 



e-x a . 



and to punch out the results as well as to store them for later use. The evaulation 
of the polynomial inTl will be faster if (2) is written in the form. 

(5) ^*(x) = |-Cn(a,+TiCa 2 +rL(a34-Ti(3 4 +na' 5 )))))^ l (x> 

To make it possible to use the LOOP order in evaluating G£ (x) this way, the 
coefficients a^ will be stored in consecutive locations in decreasing order. The 
LOOP program will be given a form applicable to an arbitrary polynomial by 
including a "dummy" coefficient a Q ss 0. Storage locations will be chosen as 
follows: 



Location 



Contents 



101-119 

200 
201 
202 
203 
204 

221 
222 
223 

301 
302 
303 
304 
305 
306 

401-410 



instructions 




(cards 1-4) 



constants 
(card 5) 



temporary 

storage 

("erasable") 



coefficients 
in <£*C*) 
(card 6) 



results 
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The program might be written as follows: 
Card Loc. Alpha. Sign A or B 



4 



n 



101 



102 



SET C 



MOVE 



005 



000 



114 



202 



C Comments 

401 Set address of 

first <£* CX) 
221 First x is x = a 



103 


NGMPY 


4- 


5 


000 


104 


EXP E 


4- 





301 


105 


MPY 


4- 


3 


000 


106 


MPY 


4- 


3 


204 


107 


ADD 


4- 


1 


000 


108 


DIV 


+ 


4 


201 


109 


MOVE 


+ 


9 


000 


110 


MPY 


+ 


3 


000 


11-1 


ADD 


— 


1 


000 


112 


LOOP B 


+ 





010 


113 


NGMPY 


4- 


5 


000 


114 


ADD 


4- 


1 


000 


115 


ADD 


4- 


1 


221 


116 


ADD C 


4 





006 


117 


TR C 


4 


8 


114 


118 


PCH 


4 





410 


119 


COND STOP 


4 





200 



000 


000 


-x^ 


000 
200 


000 
222 


e " x2 2. ^-X^ 


221 


000 


px 


201 


000 


1+ px 


000 


223 


n=r/Ci+Px) 


301 


000 


a 5 into 000 for LOOP 


223 


000 


prev. res..Tl 


302 


000 


add next coeff. 


005 


110 


loop in polynomial eval 


222 


000 


polyn. 'Cp'CX) 


201 


C 1 


(x)= 1+ prev. res. 


203 


221 


x + A =• next x 


114 


001 


next<£*(x) address 


412 


103 


test for end 


401 


410 


punch two cards 


221 


500 


end; stop, display last 
x, go to 500 on PRO- 
GRAM START 



500 Next instruction in the problem. 

An important remark should be made: If there is no shortage of storage 
locations and if the programmer does not mind writing a somewhat larger number 
of instructions, the running time for many problems can be decreased and the 
logic simplified by "unwinding" the innermost loop, i.e., by writing out the 
mathematical instructions in the loop in a straight sequence instead of using the 
LOOP operation. In the present problem, a sequence containing five MPY and four 
ADD instructions could replace the instructions 109-112 and also eliminate the use 
of the dummy coefficient a . The execution time for the polynomial loop would be 
reduced by nearly 1/3 and the LOOP operation could be used to replace the ad- 
dress change operations in the outer loop. The polynomial evaluation accounts for 
about 1/2 of the total running time of this problem. In many large problems, the 
innermost loop consumes an even larger fraction of the running time, making it 
important to program the innermost loop efficiently even at the expense of apparent 
inefficiencies elsewhere. 



The second illustrative problem reads as follows: For a given set of numbers 
v =1, 2, . . . , 50, not necessarily equally spaced, the values of the Chebyshev 



polynomials T n (x v ), n — 1, 2, 



10, are to be computed using the recursion formula 
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(6) T Tl+l Cx v )=<ax v r n Cx v )-T Tl -,Cx Y ) ; 

(T Q (x )sl, Ti (x )= x v ) and punched out in a compact form. 
In addition, the sum 

(7) fCT.oCXv)! 2 (Xy+,-Xy), 
V=lVl-X v ^ 

(Xc-i = 1) is to be punched out and the operator is to be given the option of also 
calling for the punching of partial sums of (7) at any time. 

Storage locations will be assigned as follows: 

Location Contents 



050 
051 

052 


o) 

1 ) constants (card 7) 


095-120 


instructions (cards 1-6) 


199 
200 
201 
202 


The sum (7) and its partial sums) 
T (x v ) s 1 
Ti(x v ) S x v 
T 2 (x v ) 


210 


Ti (x v ) 


300 


2x v 


301 


[VV] 2 


302 


[T 10 (x v )] 2 Vl - x/ 


400 +■ v 
451 


x v 



output 



temporary 

storage 

("erasable") 



input 

(cards 8 - 18) 



In addition, locations 1-8 will be used in connection with a trick in program- 
ming the LOOP. 

The program may be written in many ways. The following is not necessarily 
the best: 
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A number of remarks are called for, many of them of general applicability: 

(a) The C-address of instruction 101 will, during execution, run through 
the values 000-008, but the result of the instruction is always called out from 000 
on step 102. This trick makes it possible to use the LOOP BC operation instead 
of address change, which is normally required if different sets of addresses are to 
be modified during a loop. 

(b) The instruction numbering was arrived at by starting the preparation 
of the program at instruction 101 with the intention of later adding an unknown 
number of initialization steps preceding it. This speaks in favor of not starting 
a program at 001. 

(c) The stop which would normally occur each time the CONS instruction 
is reached may be bypassed when found superfluous without any sacrifice by 
turning the PROGRAMMED STOP switch to RUN, since no COND STOP, SIN or 
LOG operations (the only other ones involving a conditional stop) are used. The 
operator decision regarding punching of partial sums is made using only the sign 
switch of the STORAGE ENTRY switches. This switch does not influence LOAD- 
ING. 

(d) The quantity x is used so frequently that it was more economical to 
MOVE it into a fixed location than to apply address modification. The converse 
applies to x v + p which is used only once. 

(e) The constant 1 appears in three locations merely in order to simplify 
bookkeeping and loading, as well as changing the number of points x v in a later 



run. 



(f) An invaluable aid in determining whether the results of a calculation 
are correct is a mathematical identity which they must satisfy, and the programming 
of such checks is strongly recommended whenever it is possible. In the present 
problem, the identity 

(iii) \ "= d dx = X 



vr- 



x 



2 



is closely connected with the computation of (7) if the x v are distributed over the 
interval (-1, 1). 

(g) An alternative method of programming the outer loop, which would 
eliminate the address change operations at the expense of somewhat increased 
card preparation, would be to key punch the x v one to a card and give a READ 
order entering one x v at a time into a fixed location during the execution of the 
program. The difficulty arising from the need for x v +jon step 109 is not in- 
surmountable. 
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XI INTERNAL STRUCTURE OF THE SYSTEM 



XI. 1. DETAILED DESIGN CONSIDERATIONS 

An expert examining the program at the end of this report will ask a number 
of questions about apparent duplication, about tight optimization in one routine in 
contrast to a lack of it in another, about the choice of operations and of methods 
of implementing them, etc. This section will attempt to answer some of these 
questions and also suggest a number of changes and additions that could be con- 
sidered for a second version of the system. Additional questions and suggestions 
from readers will be genuinely appreciated. 

In the early stages of system design, the following requirements were among 
those agreed upon, in addition to the general principles discussed in Sec. I. 2: 

(a) The arithmetic operations and those logical operations most likely to 
occur in inner loops (LOOP and certain TRANSFER operations) must be as fast 
as we know how to make them, regardless of the expense in storage. 

(b) The system must occupy at most 1000 memory locations. 

(c) The special functions must have full accuracy and unlimited range and 
most of them should be as fast as these requirements and available storage permit. 

(d) Optimum programming, (see the 650 Manual of June, 1955) in addition 
to being necessary for the attainment of (a) and (c), should be used locally in any 
program where the gain is significant but not at the expense of extensive rewriting 
of previously completed programs. 

(e) The programs must be so written that if the machine stops on any pro- 
gram step in a subroutine and control is transferred elsewhere before restarting, 
the subroutine, where the stop occurred, is left in a condition which assures correct 
operation the next time that subroutine is used. This implies that if a subroutine 

is used in more than one program, it must be initialized by each program rather 
than having a normal form used in one program and temporarily being changed at 
the beginning of other programs when needed there and then restored to normal at 
the end. 

(f) To facilitate changes, the individual programs (or "decks", 1-20, see 
Sec. XI. 7) that make up the system should be as independent of one another as they 
can be without excessive waste of storage. This requirement was not fully adhered 
to near the end of the programming task. 

As a result of these requirements and of some oversights in programming, 
there are a number of storage registers which could be made available without any 
loss in system performance and a number which could be freed at some sacrifice. 
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A brief guide for finding such registers will now be given followed by a number of 
suggestions for their possible use in a revised version of the system. 

The 6 vacant 2 -code locations and the 11 vacant registers listed in deck 5 
are, of course, available. The only distinction between them is one of mnemonics 
in connection with the choice of operation codes. In addition, it appears possible 
to salvage 22 registers essentially without loss by the following substitutions, but 
a careful check followed by machine testing is advisable: 



)eck 


Card 


Loc. 


Replace by 


20 


113 


1801 


1848 


16 


48 


1240 


1138 


18 


103 


1896 


1138 


6 


70 


1360 


1160 


16 


6 


1230 


1338 


18 


77 


1887 


1137 


2 


32 


1058 


1955 


12 


60 


1639 


1289 


16 


30 


1244 


1245 


19 


103 


1702 


1103 


17 


63 


1480 


1980 


10 


36 


1331 


1358 


8 


32 


1166 


1241 


12 


56 


1423 


1674 


17 


22 


1485 


1285 


18 


97 


1842 


1504 


17 


29 


1495 


1297 


5 


11 


1252 


1952 


5 


12 


1255 


1955 


5 


13 


1260 


I960 


5 


21 


1277 


1977 


5 


23 


1283 


1983 



It is, of course, necessary to determine, by sorting on instruction and data 
address, all places where the locations listed are referred to. 

Registers that may be freed at a price in speed include, above all, nearly 
40 extra registers used in the arit%m&tic routines in calling out A and B, splitting 
them up and storing the parts. This/is done separately in each of decks 12, 13 
and 14 to accommodate minor differences that facilitate optimization. To combine 
these steps without any loss of timfe is a task which, if possible, would require 
re -optimization of a substantial pa/rt of the system. At the expense_of one revolution, 
they may be combined easily. Similarly, making the dissection of B common to 
all 0, -routines would result in a substantial saving at the expense of lost time in 
cases (such as LOOP and TREXP)| where B is irrelevant. To make this dissection 
common only to those routines where it is needed would be less profitable. 
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At some sacrifice in external characteristics, registers may, of course, be 
freed in any number of ways. If, in tracing, the modified instruction (redundant 
but convenient) is omitted, seven steps are eliminated. The MOVE operation for 
At 000 is easily programmed in terms of LOOP BC and MOVE 000 and could be 
omitted, as could the special functions in degrees and to base 10 (or radians and 
base e, respectively). 

A number of suggestions for changes and additional operations will now be 
listed. Suggestions (1) - (3) use only the vacant registers and operation codes 
listed in deck 5 and can consequently be added to the system without difficulty at 
the option of any installation or individual programmer. For temporary use, they 
may be punched on separate cards and loaded after the system deck, in the case 
of (1) and (2) and after the punching deck, in the case of (3). Such cards should 
have an x-punch in column 5 and the four-digit location in columns 6-9. 

(I) Add an 2 -operation defined as follows: 

°2 = 453 SWITCH Transfer on switch. Control is 

transferred to C if the Storage 
Entry Sign Switch is set to minus, 
to B if it is set to plus. 

This operation bypasses the stop that would occur if the same function were 
programmed by a CONS and a TR SGN order. It might be particularly useful in 
connection with tracing when it is desired to start tracing after a certain amount 
of running time has elapsed or for following the progress of a calculation by 
occasional punching of intermediate results at the discretion of the operator. 

The coding for SWITCH consists of the instruction: 



1453 10 



8000 1015 



Read console, 
routine. 



Go to TR SGN 



The execution time is 19. 2 ms. 

(2) Add an O^operation called COUNT having the same counting and testing 
properties as the LOOP orders but using a counter independent of the loop box and 
not capable of modifying instructions. Its function can be duplicated, e. g. , by a 
SET A, an ADD A and a TR A instruction. Its advantage lies in its speed and 
simplicity. The execution time is 24. ms.when control is transferred to C and 
33. 6 the last time when control proceeds ahead and the counter is reset. A formal 
definition follows: 



2 = 800 



COUNT 



The number standing in the counter 
is increased by 1. Its new value 
is compared to B. If B is greater, 
control is transferred to C. Other- 
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wise, the counter is reset to 
zero and control proceeds to 
the next instruction. The 
counter is also reset in loading. 

If COUNT is used extensively, an expansion of the tracing program to punch 
out the contents of the counter, e. g. , in place of the problem number, would seem 
desirable. 



The coding for COUNT reads as follows : 



Loo. 



Op. 



1800 


10 


1314 


10 


1323 


11 


1381 


46 


1337 


10 


1396 


21 


1391 


16 


1066 


20 


1356 


00 


1317 


00 


1194 


20 


1378 


24 



Data 

1356 
1317 
8002 
1337 
8001 
1356 
8002 
1356 

0000 
1000 

0000 
1356 



Instr. 



Remarks 




0000 
0000 

1378*1 
1178J 



Call out and increase 

the contents, N, of the counter. 

Test N+ 1- B. 

On — , store N+ 1 in the counter; go to C 
(in TR SGN program). 

On+- , reset the counter, go to General 
Int e rp r e tati on . 

The Counter 
Constant 

Change in LOADING to reset 
the counter. 



Note: If the COUNT program is loaded separately, the card loading 
zero into 1356 (step 9 in the program) must be included. 

(3) Include in the punching program (deck 9) a routine that prevents unwanted 
numbers from being punched out when the word count is less than the normal max- 
imum. This can be done on a tabulator with sufficient selector capacity (see Sec. 
XI. 5). In the 650, it requires five locations and increases the execution time of 
the PCH order by 24. ms.when the word count is less than the normal maximum. 
The program, which may be punched on three cards, reads as follows: 



1949 


44 


1306 


1095 


1306 


20 


1980 


1307 


1307 


20 


1981 


1308 


1308 


20 


1982 


1309 


1309 


20 


1983 


1044 


1044 


20 


1984 


1973 



(4) Make room for the tracing program to be on the drum in parallel with 
the punch program, replacing the mode-of-operation deck (see LOADING) by an 
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x-punch on the transfer card or a setting of the storage entry sign switch. An ex- 
pansion of the loading program (about 7 steps) or of general interpretation would 
be needed, and the present overlap between tracing and punching is 34 registers. 

(5) If MOVE is omitted, except for A= 000, make this an 2 -operation and 
use the vacant Opcode for NGDIV. Alternatively, add a fast ? -operation, "NEG", 
identical with MOVE 000 except that it changes the sign of B. If NEG were available, 
however, it might be used in cases where, by slight reprogramming, a better pro- 
gram using NGMPY could be written. 

(6) Increase the number of logical operations, adding to the flexibility of the 
system and to the confusion of the beginner: Have a register called the "address 
counter", addressed, e.g. , by 000 or by special operations and SET instructions 
referring to the address counter (as in 701 Speedcoding) where the present SET 
instructions refer to their own C-address. Have a set of TR A, TR B, TR C orders 
which automatically increase the address referred to by 1. These would have to be 
alternative to the address transfer orders in present use unless vacant 0. -codes are 
produced. • 

(7) Make use of addresses now ignored in some operations. For example, 
in CONS, use the B-address to call out a number B for console display when the 
machine stops. In TR OUT, make B a "return address" similar to that in TR SUBR. 
In START TR, or a new tracing order supplementing it, let B (or C) designate the 
number of steps to be traced before an automatic discontinuation of tracing. 

(8) Make Program Loading reset the registers below 1000 to zeros, unless 
told not to by an x-punch on the first card being loaded. 

(9) Have a conditional stop, or an operation effecting such a stop, on loss 

of accuracy in ADD and SUB, analogous to those in SIN and LOG. In many problems, 
particularly in connection with tests, such loss is legitimate, however,and a stop 
undesirable. 

(10) Replace or supplement the present error stops by the punching of an 
"error card". 

(11) Introduce an operation similar to ST TR ERAS, perhaps replacing it, 
which will cause the machine to trace the first, second and last repetition in a 
loop. 

(12) Add another LD-STD pair (at no loss in time) to General Interpretation 
(see cards 26 and 27) making ADD and SUB, as well as MPY and DIV available as 
internal subroutines. 

(13) Cut the execution time of several subroutines, such as, the arc tangent 
program, by making minor rearrangements, usually involving the expenditure of 

a few additional registers. 
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(14) Add an (^-operation, SPACE, which causes an x-punch to appear on the 
next card punched. 

(15) Interchange the functions of registers 1002 and 1702, causing the machine 
to stop sooner if a programmer accidentally attempts to continue upward from 
instruction 999. 

(16) Investigate whether a carry can ever occur on card 78, deck 18. If not, 
put the registers used on cards 78-82 to better use. 

(17) Replace or supplement the arc tan operation by an 1 -operation, ARG, 
which gives the argument (angle) of the point whose coordinates are (A, B). 

XI. 2. RELATED SYSTEMS 

Several systems supplementing the present one suggest themselves: (a) A 
system of symbolic or regional programming where the machine assigns absolute 
addresses in connection with loading; (b) A system externally identical with the 
present one, or very nearly so, operating on complex numbers, probably with real 
and imaginary parts in 8-2 floating -decimal form; (c) A system externally identical 
with the present one, or very nearly so, operating on double precision floating- 
decimal numbers, e. g. , 16-4; (d) A system of "formula translation" or "automatic 
coding" (such as, the IBM Fortran for the 704) putting on the machine as much as 
possible of the burden of translation from a set of mathematical formulas to a 
program. 

XI. 3. NUMERICAL METHODS 

The study of numerical methods for calculation of the special functions included 
in the system was not nearly as exhaustive as would have been desirable and no 
claim to an optimal choice is made. 

The square root is computed by Newton's iteration method, 

ATI 

where 1 ^ B, ^ 10, using the initial approximation 

(2) XQS1 + .22BJ . 

The evaluation of the trigonometric and exponential functions is based on RAND 
approximations (see Approximations for Digital Computers by Cecil Hastings, Jr. , 
Princeton University Press) to sin "J X and 1 X for < x < 1. Resembling the 
approximations obtainable by expansion in terms of orthogonal polynomials, these 
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approximations are in general somewhat more efficient than partial sums of Taylor 
series for a prescribed interval and accuracy, but it is not obvious that a further 
reduction of the argument followed by the use of a Taylor expansion could not have 
been better in the present case. For small x, as stated previously, the formula sin 
x = x (in radians) is used in order to retain significant figures. 

The logarithm and arc tangent are evaluated from fixed-length partial sums 
of power series after preliminary reductions of the argument, since eight-digit 
RAND approximations were not available. For the logarithm of B,, 1< B, < 10, 
the substitutions — 1 

O) U= j^ ; 7=-^ ; for B,<e, 



(4) u = fe? >V=U5 ; WB,>e, 



" I. 
e 

(5) t -^ , 

are followed by the evaluation of 



(6) log e B,= lo 5e ]i| + V = 2tD + ^ + ,,, + tl ]+V. 

The constants, V~e~and e 1 - 65 were arbitrarily chosen within the intervals that would 
lead to a minimal number of terms in (6). For x near 1, the logarithm is inherently 
less accurate than x since 



(7) A lod X=4f X 4. 



2 — X 



X 



and logxs:0 whereas x ^ 1. No substitution comparable to sin x = x can alleviate 
this difficulty. The use of a second (^-operation for log (1 4- B) was considered 
but was rejected due to space limitations. This second logarithm could not be used 
to replace the present one for all values of the argument, since if the logarithm of 
a small number, say 10 -10 is desired, the substitution, 

(8) l-i-B = lo" 10 

forced upon the programmer, yields B = -1 exactly (in the eight-digit system used) 
with all digits of the input lost. 
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For the arc tangent, the reduction (after the argument is restricted to 0< 
B ^ 1 by the obvious properties of the function) is based on the formula 

x-y 

(9) arc tan x= arc tan y +■ arc tan — -f- 

which is merely the addition theorem for the tangent rewritten. With y= .6, the 
use of (9) gives the desired accuracy in 



(10) arc tan z 






with z = ^r/ for x>.28 and z = x for x <„28. There is again some leeway in the 

choice of these constants. For small x, the substitution arc tan x = x is used to 
preserve significant digits. 

XI. 4. CONTROL PANEL WIRING FOR THE READ -PUNCH UNIT 

The control panel for the 533 Read-Punch Unit associated with the 650 is 
wired as follows: 

Col. 1, 1st Reading, to LOAD. 
R 4- Sign, jackplugged. 
P+ Sign, jackplugged. 

Col. 5, 1st Reading, to Pilot Sel. 1 X PU. 

Rd. Hold to PS1 Hold. 

Read Card C, Col. 6, to PS1 T. 

Read Impulse to PS1 N. 

PS1 C to Storage Entry C, Word 1, pos. 3 (from the left). 

Read Card C, Col. to Storage Entry C. 

Wd. 1, pos. 4-6 

Wd. 2, pos. 6 

Wd. 3, Sign 

Wd. 3, pos. 1-10 

Wd. 4, Sign 



7- 


-9 


10 




11 




12- 


-21 


22 




67- 


-76 


77- 


-79 



Wd. 8, pos. 1-10 
Wd. 9, pos. 4-6 
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Read Impulse 12 Wds. 1, 2, 9, 10, Sign 

Wd. 1, pos. 7-10 
Read Impulse ^ Wd. 2, pos. 7-10 

Wd. 9, pos. 7-10 



Word Size Emitter to Word Size Entry C 



10 



Wds. 3-8, 10 
Wd. 1 
Wd. 9 
5 Wd. 2 



8 Wd. 1 

7 Wd. 9 



Wd. 


1, 


pos. 


3-6 


Wd. 


2, 


pos. 


6 


Wd. 


3, 


Sign 




Wd. 


3, 


pos. 


1-10 


Wd. 


4, 


Sign 




Wd. 


8, 


pos. 


1-10 


Wd. 


9, 


pos. 


4-6 


Wd. 


2, 


pos. 


10 



Storage Exit C to Punch Card C, Col. 

Wd. 10, pos. 3-6 !_4 

6-9 
10 
11 

12-21 
22 

67-76 
77-79 
80 

Double Punch and Blank Column Detection as available and desired. 

XI. 5. CONTROL PANEL WIRING FOR THE TABULATOR 

The IBM accounting machine or tabulator used for printing from the cards 
associated with the interpretive system may be expected to perform some or all 
of the following tasks: 

(1) Automatic selection of different printing forms (i. e. , zero control and 
spacing between items on a line) for data cards and trace cards. 

(2) Selective spacing between lines. 

(3) Suppression of unwanted words from cards with word count less than 
5 (or 6). 

Since there are many tabulator models, it is not feasible to provide a detailed 
wiring diagram in this report. Instead, suggestions of general applicability will be 
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given. 

Exact selector requirements depend on the characteristics of each machine. 
As an example, requirements on a 416 will be given: The printing of signs requires 
6 single-position selectors with X-pickup and 6 positions of 11-12 separation (either 
special attachments or 6 positions of a selector transferred by an 11-1/2 impulse). 
Task (1) requires a digit selector (which may be put to duplicate use in task (3)) and 
34 selector positions with digit pickup (delayed pickup). Task (2) requires only 
3 one-position selectors on the 416. Task (3) requires 55 selector positions with 
delayed pickup and some comparing units or a five-position selector and a digit 
emitter for control. 

The problem in connection with task (1) is to get the desired zero control 
and spacing in the two cases with the same setting of the hammersplit levers (also 
called zero suppression levers) and hammerlocks, on machines where zero control 
is not performed on the control panel. On a tabulator with 89 type bars, this may 
be done as follows: 



Type bar 



Alpha. 



1-4 

5 
6-8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 



Direct Wiring 
Col. 1-4, II*. 
Col. 7-9, II. 



Data Cards 



Trace Cards 



Col. 5, II. 
Col. 6, II. 



Col. 11, I. 



Col. 12-21, II. 



Col. 22, I. 



Col. 23-29, II. 



Col. 11, 


I* 


Col. 12, 


II. 


Emit 10. 




Col. 13, 


II. 


Col. 14, 


II. 


Col. 15, 


II. 


Emit 10. 




CoL 1.6, 


II. 


Col. 17, 


II. 


Col. 18, 


II. 


Emit 10. 




Col. 19, 


II. 


Col. 20, 


ii. 


Col. 21, 


ii. 


Emit 10. 




Col. 23, 


ii. 


Emit 10. 




Col. 24, 


ii. 


CoL 25, 


ii. 


Col. 26, 


ii. 


Emit 10. 




Col. 27, 


ii. 
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Type bar 


Direct Wiring 


32 


( 


33 


) 


34 .. 


/ 


35 


Col. 30, II. 


36 


Col. 31, II. 


37 


Col. 32, II. 


38 




39 




40 


«•«. 


41 




42 


I 


43 


1 


44 


I 


Num. 1 


i 


2 


A 


3 


Col. 41, II. 


4 


Col. 42> II. 


5 
6 
7 


Col. 43, II. 


Col. 44, I. 


8-17 


Col. 45-54, II 


18 





19 


Col. 55, I. 


20-29 


Col. 56-65, II 


30 





31 


Col. 66, I. 


32-41 


Col. 67-76, II 


42 


. 


43-45 


Col. 77-79, II 



Data Cards 



Col. 33, I. 



Col. 34-40, II. 



Trace Cards 

Col. 28, II. 
Col. 29, II. 
Emit 10. 



Emit 10. 
Col. 35, II. 
Col. 36, II. 
Col. 37, II. 
Emit 10. 
Col. 38, II. 
Col. 39, II. 
Col. 40, II. 
Emit 10. 



(*) The symbol "II" denotes wiring from the second brushes (on some machines 
called "third reading") whereas "I" denotes wiring from the first brushes ("second 
reading") through a selector that separates ll's from 12*s to the X-PU of an X- 
distributor ("pilot selector") through the transfer point of which an emitted 10 goes 
to the type bar in question. 

The hammersplit (zero suppression) levers alpha. 4, 11, 25, 37 and num. 5, 
17, 29, 41 and the long hammerlocks alpha. 12, 16, 20, 24, 26, 30, 34, 38, 42 and 
num. 2 are raised. Left zero carry clips of width 3 are attached to hammersplit 
levers alpha. 6-8 and num. 43-45. 

Trace cards are distinguished by the presence of an 8 in column 80. This 
impulse is wired through a digit selector to the digit pickup of a row of selectors 
("class selectors" or "co-selectors with a controlling pilot selector") with a total 
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of 34 sets of points and to the hammerlock control hub. 

For performing task (2) in the manner specified in Sec-. V. 2. , four external 
switches are needed. On a 416, the single-double spacing lever and the minor, 
intermediate and major control switches can be used; on machines with four 
pluggable switches there is no problem. 

Spacing every 10 lines (alternative (b), Sec. V. 2. ) may be accomplished by 
adding 1 to a counter on card cycles and using the carry (which, on a 416, is 
automatically available at the counter total exit) to initiate a minor cycle during 
which spacing takes place and the counter is cleared. 

Spacing alternatives (c) and (d) both involve inspecting a card column at two 
reading stations and taking certain action when the digit at the first brushes is less 
than the digit at the second brushes. This may be done by wiring from second 
brushes to a comparing entry and from the corresponding comparing exit to the 
immediate pickup ( ,r ZFS P. U. u on a 416) of a selector through the transfer point 
of which the digit at first brushes is passed. In case (d), this digit is taken directly 
to cause spacing; in case (c) it is wired to the digit pickup of a selector which 
initiates spacing on the next cycle. 

Task (3) is easily accomplished if a sufficient number of selectors are avail- 
able. Since the same task can be performed on the 650 at the expense of 5 locations, 
(see Sec. XI. 1. ) the tabulator wiring will not be discussed here. 

XI. 6. SYSTEM LOADING 

The interpretive system deck, normally with 6 words to a card and an x-punch 
in 5 to get the 1 in column 6 picked up as a leading digit of the address, is loaded by 
a deck of six self-loading cards (12 in col. 1) with 7 words to a card and a card 
number in the eighth word. The first card serves the sole purpose of making a 
fixed console setting possible. The System Loading program on these cards operates 
as follows: 



8000 


70 


1951 


1333^ 


1951 


70 


0004 


0152 


0004 


70 


0053 


0152 


0053 


70 


0106 


0152 


0106 


70 


0153 


0152 1 


0153 


70 


0204 


0152^ 


0204 


70 


0251 


0056 


0056 


60 


0001 


005T] 


0055 


11 


0251 


0155 


0155 


44 


0152 


0102J 


0102 


60 


0002 


000"7^ 



Read in the loading pro- 
gram from six load 
cards (B into 0001, C+ 1 
into 0002). 

Read a system card (non-load). 

Go to stop if expected loc. 
B'^t loc. on card, L. 



■61- 



0007 
0205 
0057 
0051 
0052 
0103 
0156 
0206 
0054 
0203 
8002 
8003 
0151 
0201 
0207 
0101 
0105 
0107 
0104 
0157 
0005 
0003 
0152 

0001 
0002 
0006 
0202 
0154 



11 
11 
46 
21 
60 
10 
15 
10 
21 
11 
69 
24 
15 
10 
11 
44 
10 
60 
44 
60 
10 
21 
69 

00 
00 
24 
69 
00 



0001 
0252 
0152 
0080 
0006 
0251 
0202 
0252 
0070 
0252 
C02531 

Cl] 

0154 
8001 
0070 
0105 
0070 
0080 
0157 
0001 
0252 
0001 
7777 

B 
C 1 
0000 
0253 
0001 



oooo\ 

0000 J 
0151 V 
8003 I 

oooqy 



Go to stop if C+ 1 <, 
B*+ n. 



Prepare accumulator for 
move. 

Store test constant. 



Move one word. 

Increase addresses by 1. 

Test for end of moving. 

Return to move another word. 
If C + 1 = B 4- n, end of load- 
ing; go to console. 
Increase location by word 
count, go to read next card. 

Error stop. 



Constants 



A similar program is used for punching out the system in condensed form in 
case extensive changes, entered on self-loading, single-instruction cards, have 
been made. 

The Reset Memory Card, mentioned in Sec. V. 3. , is a load card with eight 
words. The program, essentially identical with one supplied by the IBM 650 Sales 
Research Group at Endicott, runs as follows: 



8000 


70 


1951 


1333 


1951 


69 


8000 


1953 


1953 


24 


0000 


1954 


1954 


69 


1957 


1955 


1955 


24 


0999 


1956 


1956 


61 


1958 


8003 


8003 


20 


Qoooi] 


0999 
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(ll 1952 8003 

0999 (00 0000 0000 

0000 70 1951 1333 

1952 00 0001 0000 

1957 11 1952 8003 

1958 20 0001 0999 

XI. 7. PROGRAMS 

The complete programs of the system are listed on the next 21 pages, (i) - 

(xxi). In many cases, but not always, a constant used in two programs is listed 
in both. 
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650 INTERPRETIVE SYSTEM. 



DECK CARD LOC. 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 

1 
1 
1 

1 
1 

1 
1 



1 

1 



2 
2 
2 
2 

2 
2 
2 
2 

2 
2 
2 
2 



1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 



16 
17 

1 8 
1 9 
20 

21 
22 
23 
24 



1 

2 
3 
4 

5 
6 
7 
8 

9 
10 
11 
12 



1095 
1014 
1043 
80 01 
1107 

061 

01 

02 

03 

04 

04 

10 

06 

6 



9 
6 
3 
1 
8 
5 
3 
7 



03 



1112 
10 46 

1076 
1034 
8 003 

1098 
1024 
1020 
1037 



2 5 1017 

26 1729 
2 7 1735 



OR 

60 
10 
21 
60 
4 6 
30 
20 
60 
30 
20 
60 
30 
44 
10 
69 

11 
61 

60 
10 
65 

60 

69 
65 



69 
24 



DATA 

109 8 
102 4 
109 8 

[ 1+1 ] 
1112 
00 
102 
8 00 
000 
104 
8 00 
000 
106 
102 
8 00 



INSTR. 



1. GENERAL INTERPRETATION. 



3 
3 
3 
3 
5 
3 
3 
7 

2 



1 1 4-) 

17 2 9 [-Increase i to i+1. 
8 lj (see cards 26 & 27 below) 
110 7 Start dissection system instr. 
Shall the loop box -be used? 



106 
101 
02 
03 
04 
04 
10 
06 
07 
800 
108 



1 

9 
6 

3 
1 
8 
5 
3 
6 
3 
1 



1017 
800 3 

8 00 1 

103 7 

104 5 

0001 
8 00 2 
104 5 



104 6 
1061 



Store C. 

Store B. 

Is an 02-operation called for? 

Go to an Of -subroutine « 

Add the loop box to the instr. 
Remove the minus sign. 



103 4 
800 3 
[1000+02 ] Go to an 02-subroutine< 

110 7] 

V Constants 
108 



08 1 
00 OJ 



1732 1735 
153 9 104 3 



28 17 32 35 0002 1445 



100 
1119 
8002 
1096 

1200 
10 59 
8 2 
1016 

1201 
1015 
1120 
10 69 



30 
1 5 
60 
69 

30 
15 
60 
01 

10 
4 6 
60 
60 



000 3 
10 30 

1 B I 
9 99 9 

3 
116 2 
[ B ] 



10 9 
1120 
102 3 
8 00 2 



Loop box (initially 0) 

Restore the multiplication 
routine (0j-3) to normal 
(The trig, program makes 
special use of it). 

.2. TRANSFER OPERATIONS . 

11191 0g»000 . UNC STOP 

8 2 - Call outTTfor console display 

109 6J 

112 Stop. 

02,-200 . COND STOP 

Call out B for console display 

Stop| Go to C on PROGRAM START 

] °2 = 201 . TR SGN 

j Test the sign of prey, result 



10 5 9 
80 2 
1016 
112 

1015 
106 9 
102 7 
102 7 



Go to C 
Go to B 



if -. 
if +. 



13 1202 60 1009 1013 



02-202 . TR EXP 
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2 
2 
2 
2 
2 

2 
2 
2 

2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 

2 
2 
2 



2 
2 
2 
2 
2 



3 
3 
3 
3 
3 



14 
15 
16 
17 
18 

19 
20 
21 



2 
3 
4 
5 
6 
7 
8 
9 

1 
2 



33 
34 
35 



37 
38 
39 
40 
41 



1 
2 
3 
4 
5 



1013 
1021 
1 029 
1036 
1012 

120 3 
1027 
1035 



1 

1 
1 
1 
1 
8 
1 
1 
1 
1 
1 



04 
27 
40 
51 
61 
02 
75 
33 
9 2 
49 
58 



1205 
1128 
8003 



36 1454 



1 

1 
1 
1 
1 



114 
145 
038 
162 
030 



1 01 
1064 
1071 
10 77 

1135 



3 
6 
3 

1 
4 

6 
3 
1 

6 
3 
6 
2 
1 
6 
3 
6 
1 
3 
2 




7 

6 
6 





5 

9 
2 
5 
7 


5 
5 
1 



65 
3 5 



61 
21 
60 
60 
60 



30 
69 
30 
10 
10 



0002 
8002 
0001 
104 5 
109 5 

1023 
0003 
1038 



102 
000 
114 
105 
111 
[ c 
000 
800 

104 
000 

[ c 



102 3 
000 3 



000 5 
1017 
0002 
8 001 
113 8 



} 



Get the exp. of prev. result 



102 1 

102 9 

103 6 
1012 
1120 

102 7 ___ __ 

10 3 5V Replace i+1 by C in 

1 4 3 J general interpretation. 



Compare it to Bj Go ahead 
if exp. < B, to C if exp. £ 



B. 



'V Rei 



■203 . TR 







5 
6 



112 

104 

11 

11 

800 

107 

113 

109 2 

104 9 

105 8 
12 3, 



1 
1 
2 
5 
3 r 



2 -20/t . TR SUBR 



Set the C-address of the 
V instruction at C equal to B; 
30 to TR to C. 



112 81 
800 3 
C C ] 

1095 



02-205 . TR OUT 
• Go to C in machine language 



107 5"^ 



120 
110 
101 



109 6) 



Og-454 . NOOP 

Go to general interpretation. 
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3. LOOP OPERATIONS 



1064) 02-001 . LOOP C 
10 7 1 Move B to the C-address 
1 7 7 > position; Call out and 
113 51 increase the loop box. 
109 3J 



3 
3 
3 
3 
3 
3 

3 
3 
3 

3 
3 
3 

3 
3 



6 
7 
8 
9 
10 
11 



1 
1 

1 

1 

1 

1 



2 
3 
4 

5 
6 
7 



1 
1 
1 
1 
1 
1 



093 
051 
54 
062 
155 
113 



18 
1 9 



1100 
10 57 
10 2 5 

1101 
1008 
1074 

1110 
1121 



11 
46 
10 
21 
16 
20 

3 
10 
10 

30 
10 
10 

69 
30 



80 
10 
80 
10 
80 
10 



0001 

1017 
103 1 

001 
1017 
1031 

1017 
0001 



I 5 11 Common LOOP steps 

115 5?- Compare the loop box to Bj 
10 6 2 if less, loop again (go to C), 
112 0j 

I I 1 31 If equal, reset ihe loop box, 
1 9 5 J go to general interpretation 



1057 
10 2 5 
1 9 3 

100 8 
1074 
113 5 

1121 
102 8 



Qg-100 . LOOP A 



02-101 . LOOP AC 



Qg-110 . LOOP AB 
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3 
3 

3 
3 
3 
3 

3 
3 
3 

3 
3 

3 
3 
3 
3 
3 



4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 

4 
4 
4 
4 

4 

4 
4 
4 
4 
4 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 



2 
2 

2 
2 
2 
2 




1 

2 
3 
4 
5 



26 
27 
28 



2 
3 

3 
3 
3 
3 
3 



9 


1 
2 
3 
4 
5 



1 
2 
3 
4 
5 
6 
7 
8 
9 

1 



12 
13 
14 
1 5 

1 6 



1 
1 
1 
2 
2 
2 



10 28 
1136 

1111 
1171 
1078 
1185 

1010 

13 4 3 
1003 

1011 

118 

1138 
1031 
10 47 
1039 
1139 



10 8001 1136 
10 1039 109 3 



1153 
1312 
1322 
1330 
1388 



1 
1 
1 
1 
1 
1 



367 
328 
340 
363 
172 
336 



69 1017 

30 0001 

10 8001 

10 113 9 

3 4 

10 104 7 

10 1017 



30 
10 



6 
4 
3 
1 
2 
3 
1 
4 
1 
1 
4 



10 
30 
10 
16 
15 

10 
66 
11 
30 

11 
16 



000 4 
113 8 



0100 

0100 
0100 



8 
1 

1 
1 

1 
1 
8 
1 
1 



8002 
000 3 
132 5 
128 3 
114 2 



102 
8 
104 
000 
132 
114 




117 1 
107 8 
1185 
109 3 

134 3 
100 3 
109 3 

118 
134 3 

00 1^\ 

100 
100 
100 1 



100 
118 
171 
132 
118 
109 
132 
113 
119 
1 09 



\ 



02-111 . LOOP ABC 



02»010 . LOOP B 



02=011 . LOOP BC 
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4. MOVE 



115 3 

13121 

132 2 

133 
138 8 
119 7 



132 
134 

136 3 

117 2 

133 6 

119 7 



Ov 



Test 



% MOVE 

N (i.e., A); go to special 
move if N-0 



Store N-l 



Go to stop if N-l+C ^ 1000. 

Test for upward or downward move. 



<■ Initialize for downward move. 



8) 




I 



Initialize for upward move. 



6 
3 

l\ 
l\ 

r\ 
9 ■ 
3J 

1} 



Store test constant * (see next p.) 
Complete initialization 



117 

800 

800 

1 1 2J- Move one word 

115 

109 

132 

117 

800 3j 

105 

136 



Test for completion 



Change addresses by 1, 
return to move next word 
If N-0, Get B and go 
to storing routine. 
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4 
4 
4 

4 
4 
4 



5 
5 
5 
5 
5 
5 

5 
5 
5 

5 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 



6 
6 



34 
35 
36 

37 
3 8 
3 9 



1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 



2 
3 
4 
5 
6 

7 

8 
9 


1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 



31 
32 



3 
3 
3 
3 
3 



1362 
8 02 
1315 

1325 
1142 
1065 



1 5 
60 
69 

69 
24 
60 



106 5 

[ B ] 
2 22 2 



80 21 
144 5J 



1 1306 77 7 77 7 



1 
1 
1 
1 

1 



307 
308 
309 
4 53 
800 



1009 
1023 
1045 



1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 




2 
5 

4 
5 
7 
8 

2 
4 
7 
8 
3 
5 
9 
1 
3 
4 
7 



77 
77 
77 
77 
77 

55 
55 
55 

88 
88 
8 8 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
8 8 



777 
7 77 
777 
777 
777 



5 55 5 
555 5 
5 55 5 

8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
888 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
8 88 8 
888 8 
888 8 
8 88 8 



10 9 5 Error stop for MOVE, READ & PCH. 



800 2 
110 2 
144 5 



} 
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ERASABLE AND VACANT LOCATIONS. 



Vacant 02-code3 



55 5 5] 

5 5 5 5 I- Inter-subroutine storage 

555 5) 

88 8 81 

8888 1 

883 8 

88 8 8 

88 8 8 

88 8 8 

88 8 8 

888 8 

888 8 

888 8 

888 8 

888 8 

888 8 

888 8 

88 8 8 

888 8 

8888 

88 8 8 

8 8 8 8 V Erasable 

88 8 8j 



INITIALIZATION CONSTANTS 
FOR MOVE 




Downward 


Upward 


■* 


B+N-l 


-B 


#* 


B 


-(B+N-l) 


### 


C 


-(C+N-l) 



30 1356 77 7 777 77771 



314 
31 7 
323 
378 
381 
337 
391 



77 
77 
77 
77 
77 
77 
77 



77 
77 
77 
77 
77 
77 
77 



777 
77 7 
777 
777 
777 
777 
777 



V Vacant registers 



38 1396 77 77 7 7 7777 

39 1066 77 77777777 

40 1044 77 7777 1111) 



1 1087 69 114 11 



10 88 69 1341 1195 




ADDRESS CHANGE OPERATIONS. 



'I* I& 5 



-v- 



3 1089 69 1342 1195 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 
6 
6 
6 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 
6 
6 
6 
6 



6 



4 
5 
6 

7 
8 
9 


1 
2 
3 
4 
5 



16 
17 
18 
19 

20 

21 

22 



2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
4 
4 
4 
4 



4 4 



1 
1 
1 
1 
8 
1 
1 
1 
1 
1 
1 
1 



19 5 
53 
11 



69 
02 
98 
59 



1140 
1341 
1342 
1072 

100 5 

10 50 



1 
1 
1 
1 
1 
1 
1 
1 
1 



4 
5 
6 
7 



1 

2 

1 

1 

2 

3 

1 

1 

1 

80 

128 

1 

1 



1 

1 

1 

1 

1 

1 

2 

2 



63 
69 
77 
8 4 
92 



9 
7 
4 

2 
9 
08 
22 



3 
4 
5 
1 
1 
3 
9 
7 
6 



1 
1 
1 
1 
1 



103 
116 
000 
160 
104 



50 1006 



51 1007 



2 
6 
3 
1 
6 
3 
6 
3 
6 
3 
1 
4 



8 
2 
3 
2 
1 



119 
8 00 
000 
107 
[ A 
[1,4,7] 
8002 
000 
8 00 
000 
104 
120 



0i»8 . TO C 

Common steps . 

Set amount of shift for TR A, B or C. 



Get the instruction 
located at A. 



a) 



3 
3 
3 
5 
3 



105 3 

1311 

1369 

800 2 

119 8 

135 9" 

11171 Separate out its A- f B-, 

112 5 or C-address. 

118 3. 
139 2 

119 9 
109 5 



Compare this to B. 

Go on if equal, to C if unequal. 



35 

35 
35 
67 



0001 
000 4 
007 




Constants 



65 8003 1163 Og-005 . SET C 
65 1103 1163 Q 2 «050 . SET B 



1500 65 1104 116 3 P «500 . SET A 



9 
2 
9 
2 
5 

9 
2 
5 






6 
6 
5 
5 
5 
5 
1 



116 

27 4 
08 

28 9 
04 
00 
16 
26 
07 

[0,3,6] 

126 5 
8 
00 
8 00 
111 
102 
102 
00 
126 
[0,3,6] 
[ B ] 



5 
3 

7 
3 

] 



3 
3 
3 
5 
3 
3 
3 
5 



Common steps 

Set amounts of shift. 



Set the address for storing 
the modified instr, at B. 
Get the instruction (!) 
located at B. 

Store the right end of B. 

Destroy the old A-, B- or 
C-address 



Replace it by C 



Attach the right end; 
store at B. 



3 5 
30 

21 



65 
66 



000 3 

000 6 

102 3 

102 3 



^ 



1267 
110 8 
109 5 



- Constants 



107 9 02-006 . ADD C 
10 7 9 Op-007 . SUB C 



-vx- 



6 

6 

6 

6 

6 

6 

6 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 

6 
6 



7 
7 
7 
7 
7 
7 
7 
7 

7 
7 
7 
7 
7 
7 
7 
7 



8 
8 

8 
8 



5 2 1079 30 0007 1004 Common to both 

53 1060 65 1023 1129 0g=060 . ADD B 

54 1070 6 6 1023 1129 02=070 . SUB B 

55 1129 30 0004 1004 Common to both 

56 1600 65 1023 1189 0?"600 . ADD A 



57 
58 



5 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 



9 

1 
2 
3 
4 
5 
6 
7 
8 
9 



1 
2 
3 
4 
5 
6 
7 
8 

9 
10 
11 
12 

13 
1 4 



1 

2 

3 
4 



1700 
1189 



6 6 10 2 3 
3 001 



04 
68 
99 
57 
13 



20 
65 
30 
69 



8 
1 3 
1 1 
1156 
1267 



2 
1 



60 
46 




1 



70 1360 
7 1 1123 
7 2 1073 



80 
1 3 



1 
1 

1 
1 
1 
1 

1 
1 
1 
1 

1 
1 



00 
33 
2 4 
94 
78 
75 
44 
4 2 



377 
376 
397 
345 

18 
12 6 



21 

21 
60 
60 



70 
60 
24 
20 
24 
24 
60 
4 5 

69 
24 
69 
24 
6 
1 



1 
1 

[ 



260 

04 5 

00 

36 

26 

12 

B 

15 
260 
26 

B ] 



3 


7 

3 



1156 



1951 
1242 
127 8 

1017 
198 6 
1952 
136 4 

195 9 
198 5 
197 6 
1061 
1951 
10 38 



15 1038 60 

17 1242 70 1951 



118 9 02°700 . SUB A 
10 4 Common to both 



106 
139 
13 5 
131 
137 
800 2] 
134 8 J 
110 6" 
126 7 
126 7 
109 5, 

109 5' 
134 8 
12 8 9. 



Common to all Agp & SUB ap. 
Store shifted ±Cj 
Set address for storing 
modified instr. at B 

Get B 



Add ±C to j B| j store 
(with original sign) at B. 



Constants 



;} 



;} 



7. PROGRAM LOADING . 
Read the first program card. 
Set an instruction for return 
from the READ routine 
Reset 000 

Reset the Loop box 
Reset the card counter 
Test the word count, n. If ^ 0, 
7 J go to READ to move n words 

If n=0, put the problem no. 
into the punch band, 
decide if tracing is required, 
(dep, on contents of 1976*), 
store the address of the 
first programmer instruction, 
go to execute it (gen. int. ) 

|30 0003 1019 if not tracing 
*121 1980 1386 if tracing 
Constant (return from READ) 



1401 
1131 

14 00 
1174 



01 
60 

69 

24 



1318 
12 7 8 



113 1 
14 4 

117 
14 8 6 



il 



8. READING OPERATIONS 
Qg-401 . READ CONS 



8000 144 5} Stop. On start, read console 
117 4") 



switches, go to storing routine. 

2 =^00 . READ 

Set return instructions, go to 



-vii- 



8 

8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 
8 

8 
8 
8 
8 



9 
9 
9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



3 3 
3 4 
3 5 
36 



5 1486 69 1241 116 8j Steps common with PCH 



Return from common steps. 
Read a card. 

Go to stop if loc. L on the card 
differs from progr. first loc» 
(the stop is in MOVE) 
Go to stop if B'+n > C+l 
(n=word count, B' "current 
first location) 



Calculate and store 
test constant 

Complete initialization 

Move one word from the read 
band into place 
Increase addresses by 1. 

Test for end of movingo 

Return to move next word 

If B'+n - C+l, end of READ 

Add n to B' to get new B 1 , 
go to read next cardo 

Constants 



9. PUNCHING . 
02*410 . PCH 




2 
3 

4 
5 
6 
7 

a 

9 

o 
l 

2 
3 
4 
5 
6 
7 
8 
9 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
2 



1188 
1318 
1241 
1324 



69 
60 
21 
24 



195 3 
128 5 
127 7 



80 31 
139 
173 6 
132 OJ 



1 1410 69 196 4 1168 



1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



168 
081 
232 
36 1 
319 
134 



24 
65 
10 
30 
1 5 
20 



72 
89 
98 
66 
73 
88 
72 
62 
99 
61 
68 
48 
7 



2 
1 
4 
6 
2 
6 
1 
1 
4 
6 
3 
1 
1 



1272 
102 3 
104 5 
3 
102 4 
129 1 

197 7 
8 00 2 
196 6 



10 8 1^ Common with READ 

12 3 2 

1 3 6 1 1 Prepare for testing 

1 3 1 9 f and initialization. 

113 4 

127 2 J 



96 9 

97 8 

29 1 
97 7 



197 8 

197 5 

8 1 

4 

196 3 

197 4 



} 



98 
9 9 
31 
97 
98 8 
9 7 21 
9 6 2 
99 
96 
96 
94 8 
97 
9 5 OJ 



9 
1 
81 



FCH program only 

Go to stop (in MOVE) if 

B > C+l. 

Set the word count 

Let B' "first loc. not yet 

punched. If B'+5 < C+l, 

go to punch 5 to a"" card. 

If not, go to further testing. 

Prepare to move n words 
to the punch band 



-vlii- 



9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 



10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 










2 1 



1 

1 

1 

1 

1 

1 

1 

10 

10 












2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

4 
4 
4 
4 
4 
4 



3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 

9 

1 
2 
3 
4 
5 



1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 8 

29 

30 



1950 
19 87 
8 002 
8 003 
199 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 



1 
2 
3 
4 
5 
6 
1 
7 
7 

7 
6 
5 
5 
9 

6 

4 
9 

3 
4 
8 



6 

8 

4 

4 

4 

4 

50 

61 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

8 002 

1968 

196 

1 97 

197 

8 

197 



1 
1 
1 
1 
1 
1 
1 



88 
89 
90 



1 
1 
6 
2 
6 
1 
1 
1 
1 
1 
1 
2 
7 
6 
1 
2 
1 
4 



1 
1 



ri 97 9 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
8 
1 



97 7 

990 

B' 

97 

98 

02 

2 

02 

2 

2 

02 

98 

97 

97 

97 

97 

00 

97 



30 000 3 

21 197 7 

5 

24 197 8 

69 

24 1978 




8 02 
1991 
19 9 2 




198 7 
800 2 

8 3} Move one word 

[1990+nlj 

Go to punch 

Return to move 
another word 

Increase the card no. 

Punch a card 

Set the location no. 

for the next card 

If B'«C+1, punching 

is completed. 

[This const, keeps trace orders 
{inoperative during non-tracing. 
(This const, will be 00 0000 1095 
[during tracing (PCH Inoperative) 
99 9 9y Constants 

199 1/ 

8 3 (When the PCH progr. is loaded 
19 9 lj- with program loading, this instr. 
.will load into 1963 (!). 

10. TRACING 

Expansion of gen , int. 

Store modified instr. for trac. 

Store C 



Store B 

Call out A if Oi / 

Store A or for tracing 

Store B for tracing 

Store prev. result for tracing 
Store loop box for tracing 
Store instr. no. for tracing 

Store original instr. for tracing 



-ix- 



1 
1 
1 
1 

1 
1 














10 
10 

10 
10 
10 
10 
10 
10 

10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 



1 
1 
1 
1 
1 
1 
1 
1 



11 



12 
12 
12 
12 
12 



1 
1 
1 
1 
1 

1 
1 
1 



3 
3 
3 
3 

3 
3 



1 
2 
3 

4 

5 
6 



37 
3 8 



3 
4 
4 
4 
4 
4 



4 5 
4 6 
4 7 

48 



4 
5 
5 
5 
5 
5 
5 



1 
2 

3 
4 
5 
6 
7 
8 



1 

2 
3 
4 



6 
7 
8 
9 
10 



1993 
1994 
1995 
1996 

14 50 
1331 

14 51 
1358 



1 
1 
1 
1 
1 



452 
332 
091 
165 
334 



8 02 



1 

1 

a 
l 
l 
l 
l 
i 
i 
l 
l 



5 8 



1 
1 
1 

1 
1 
1 
1 



4 45 
404 
414 
427 
430 
4 39 
4 47 
8 02 



1082 
1492 
14 49 
14 57 



15 
20 
71 
60 

69 
24 

69 
24 

69 
24 
65 
16 
69 
24 

09 
60 
6 9 
04 
36 
30 



21 
44 



2 
2 
6 
4 
3 
1 
6 
2 



9 1443 24 



24 
65 
30 
1 5 



1024 
198 6 
197 7 
126 

1976 
1061 

115 4 
1061 

1976 
1061 
109 8 
1173 
1141 
[i + l] 



5 400 

000 3 
00 6 

198 
001, 
196 5 



19 9 4V Store card no. for trace cards 



100 
102 
143 
000 
144 



C c ] 



128 9 
104 5 
00 3 
1460 



10 83 
1542 
1599 
120 7 
8 02 



11 14 7 3 

12 148 3 

13 1493 



24 
65 
30 
1 5 
61 

21 
65 
4 5 



128 9 
104 5 
3 
1310 
C B ] 

127 8 

128 9 
149 6 



199 5. 
199 6 
110 5 



Punch a trace card 
Return to normal gen. int. 



13 3 1) 02=450 . START TO 

1 9 5 f Modify gen. int. to include 

tracing expansion 
13 5 81 02=451 . STOP TO 
109 5; Restore gen. int. to normal 

13 3 2] 02-452 . ST TO ERAS 

10 9 1} Modify gen. int. to include 



116 51 
1334 
800 2 
109 5j 



tracing expansion. 

Replace the present progr. 
instr. by a NOOP 



88 
96 




1 
197 2 



41 
8 




1 

1 
1 



196 8 




> Constants 



Word count ^ trace identif . 
Make PCH inoperative when tracing. 
Will be 30 0003 1019 when not trac. 
When the trace program is 
loaded with program loading, 
this will load into 1963 (!). 
11. STORING THE RESULT 
(Common to all math, routines ) 
Store C in 000 for prev. res. and 
in 1009 for cond» transfer ops. 
If C-000, go directly to gen. int. 

If C ft 000, store C in C, 
go to gen. int. 



10 9 5 Constant 



149 2 
144 9" 
14 5 7 
800 2 



12. ADDITION AND SUBTRACTION 
0t«1 . ADD 

Get B 



5 8002 60 [ B] 147 3. 



15 4 2 0i:"2 . SUB 
159 9) 



120 7 
800 2 

147 3. 

148 3 

149 3 
149 7 



y Get -B 

Common steps 

If A-000, get -R directly. 



-X- 



1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 



Get A if A / 000 




Split up 
Split up 



Put |b-a| into shift instr. 



8 V If a-b, form Ai+B-|»CJ, 

9 

5J 

OVlf |C|| > 10, C] = CJ/lO. 

ll If Ci < 0, get c«a+l, 

1 r combine with Cj , 

7 ) go to store. 

5. 

l\ Same if C] > 

1] 

5 If Ci»0, go to store machine 0. 

If [C{ | < 10, get Ct by 
shifting; correct exp. 

3 If a / b, prepare to shift. 
5 

3 I If_|b-a| > 10, get C-A or 

5 f OB, go to store, 

5J 

3 

2" 

- If b < a, shift Bf . 



If b > a, interchange 

Ai and B-\ t replace a by b. 



12 
12 
12 
12 



6 6 
6 7 
6 8 
6 9 
70 



1310 
14 6 
140 9 
1433 
1456 



61 
6 
60 
31 
01 



147 3 
147 3 
144 6 
14 13 



Constants 



-xi- 



13 
13 
13 
13 
13 

13 
13 
13 
13 
13 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 

3 

3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 



1 
2 
3 
4 
5 

6 
7 
8 
9 
10 



1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 

3 
3 

3 
3 
3 
3 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
5 
5 
5 
5 



1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 

6 

7 

8 

9 



1. 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 



10 84 
1592 
1549 
1557 
8 02 

10 86 
1692 
1649 
1657 
80 02 



1 
1 
1 
1 
1 



80 
1 5 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



8 
6 
3 
6 
5 
2 
7 
5 
1 
5 
3 

7 
3 
7 
6 
5 
3 
1 



88 
91 
46 
64 
48 
59 

1 7 
73 
77 
33 
28 
84 



24 
65 
30 
1 5 
60 

24 
65 
30 
15 
61 



28 
81 
32 
39 



2 
6 
4 
3 
1 
6 
6 
3 
4 
2 
2 
6 
3 
4 
2 
6 
1 
1 
2 
4 
4 
6 
1 
3 
6 

3 

4 

6 

3 

6 

3 

3 

4 

4 

1 

1 

3 

6 

30 

4 6 

18 

1 7 

3 5 



1 
5 
5 

5 



4 
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6 
4 
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128 9 
104 5 
3 
1560 
[ B ] 

128 9 
104 5 
3 
1660 
C B. ] 



28 3 
28 9 
59 6 
00 3 
50 8 
A ] 



00 

1 51 
2 
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5 
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8 00 
156 
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8 00 3 

1 
167 7 
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00 
128 
000 
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125 
150 
8 00 
8 
00 
8 
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128 
128 
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155 
800 
157 



169 2 

164 9 

165 7 
800 2 
157 8 



13. MULTIPLICATION 
Oi-3 . MPY 

V Get B 



Oi»5 . NGMPT 
' Get -B 

Common steps 
Get A 

If AtO, C-0 directly. 
A»Ai,a. Split up and store 
B»Bl,b. Split up and store. 



Calculate and store 
49 - (a+b) 

Go to error stop if exp. of 
prod, would be out of range 

CJ - At • Bf 



Common with DIV. 



Cj/10. 



144 5. 



If Cj > 10, C 7 

If CJ < 10, prepare to 
correct exponent 



Jlf Ci«10 (due to carry in 
[rounding), go to correct. 



Calculate c 
attach it to 



* a+b-49-corr. j 
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13 54 1595 6 9 5555 1338 Error stop, for MPY and DIV. 
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1 338 
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8 
64 



3 3 1610 
3 4 1659 
3 5 1654 



1673 
1527 
1374 



1350 
1169 
8002 
1137 
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49 
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30 
15 
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67 
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8 03 
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04 
00 
61 

B 
28 
28 
64 
00 
65 

A 
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26 



28 
29 
26 



128 3 
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9 

] 

2 
7 
4 
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2 
7 
3 
2 
8 
4 
1 
1 
7 
4 
7 
4 
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1 
8 



60 
60 
50 

4 5 1626 
4 5 168 
69 3 33 3 



000 3 
1022 
[ B ] 
8 002 



157 8 
157 8 
165 5 



V Constants 



B=0 tested below 



80 
16 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



00 1. 



14 4 5 Go to store after error stop 

15 3 9 
U. DIVISION . 

7^-Get B 
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31 
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6 
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3 
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3 
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3 
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68 5 

64 4 

5 1 
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2 1 
9 5 
7 8 

3 8. 

4 5 
1 7 



VGet A 



• A-Ai,a. Split up and store 
Aj-6 tested below 



- B"Bi ,b. Split up and store 



Store shifted Aj, 



* Calculate and store b-a-50. 

|Go to stop (in MPI) if exp. of 
J quotient would be out of range, 

Ci=Ai/Bij go to exp. corr. 
(in MFY) and storing routines. 



16 7 31 
160 6 



- Constants 



13 7 4 Go to stop if B-0. 

14 4 5 set C-0 directly if At 
13 3 8 Error stop for B-0- 



1169^ 
800 2 
113 7 
144 5J 



15. ABSOLUTE VALUE . 
Og-350 . ABS 

Get |b| - C, go to store it. 



-xiii- 



15 



5 1022 6 7 



113 7 Constant 
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16. SQUARE ROOT 
02=300. SORT 

B 

Stop if B < 0. 



If B-0, 00 directly. 
B-Bi,b. Split up and store 



x =l+.22Bi 



Calc. and store Bi/x^. 



If iBi/xj, - xj < 10~ 8 , go to end 



x^ - 1/2 to/xn + x n ) 



- Calc. 1/2 (b+50). 



b is even, go to end. 



b is odd, Ci » Xjj -/LQ. 



Equip with exp. , go to store 



b is even, Ci-Xjj. 



1 
1 
1 
1 
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419 
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Constants 



51 1150 03 1622 776 6J 
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17, EXPONENTIAL 

1 1301 6 9 1514 1467 ? °301 . EXP E 

2 1351 69 1464 1467 Op°351 . EXP 10 
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5 58 
61 7 
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40 
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602 
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520 
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00 
148 
44 
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00 
26 
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52 
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[ |b-49 1] 
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27 8 
58 7 
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27 



3 
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8 00 
1 50 

127 

800 3 

165 8 

127 



8 00 
161 
127 
8 00 
156 
1270 
8 00 3 
151 
1 27 
8 00 
146 3 
1270 
8 00 3 
166 5 
8 3 



3 
2 

3 
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3 
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142 5^ 
148 8 
800 2 
1553 



Common steps . 

Set an instr. for EXP E or EXP 10. 

Get B°Bi ,b. Split up 
and store. 



Set an instr. distinguishing 
between pos. and neg. exp 



Form b-49. Go to special 
routines if |b-49| £ 10. 



Set instr. for left or right 
shift of |b-49|. 
Go to EXP E routine 
Go to EXP 10 routine 

For eB-lOwB, calc. |mBi|. 

For 10 5 , get |B||. 

Shift |S|-n+\. (0£\< 1). 
Store \. 
Store n. 

If |B| > 50, go to stop. 



lO^ - [l+a 1 \+.,.+a ? \ 7 ]2 
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8 00 3 
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For 1 > 0, round lo\ 
equip with exponent n+50, 
go to store. 



r\ 



For B < 0, calculate 10 



Equip 10 with exponent 

49-n (if 10"^ < 1) or 

50 (if lO-'X. - 1); go to store. 



Go to 8 top if b > 59, 
set Ol if b £ 39. 



155 1 
150 2 
155 3 
168 7 
143 8 
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48 

00 

26 
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173 

35 7 

08 8 

27 7 

005 
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y Constants 



9 6 1393 69 6 66 6 1338 



Error stop for exp. , log. & trig, 
18. LOGARITHM. 



1 1302 69 1807 1818 Op-302. LOG E 

2 1352 69 1815 1818 Qg-352 . LOG 10 

Common steps ,. 



Get B-Bi ,b. Split up, 
store exponent 
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1 39 3". 
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139 3 


175 1. 
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17 51 


60 
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18 9] 
18171 
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1 2 


18 9 


11 


1812 



Go to stop (in EXP) if B 7 ^ 0. 
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18 2 1 8 7 1 y If B, < e, set u - B^-Ze, 

v-,5. 
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If Bi > e, set u » Bi/e 1 * 65 , 
v=1.65 



u-1 
u+1 
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j-1-i^.g* 
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182 6 
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18 2 2 
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End of common steps 



For LOG E, store L, 
calculate Mbj. (by 



■b-50) 



8 4 6] 

810 
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8 8 

89 3 

84 

44 

86 

86 

87 



For LOG 10, store mL, 
calculate b 1# 



Common steps 



■(Go to stop on loss of more than 
(.two digits of accuracy 
jGo to store machine zero 
(.if C has seven zeros. 
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8 75 
834 
8 43 
8 04 
811 
8 78 
8 08 
8 79 
8 87 



65 
31 



1750 
18 54 
1874 
1872 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



8 
81 



7 
5 

836 
812 
8 28 
867 
8 63 
816 
8 30 
8 86 
869 
832 
848 
806 
8 42 
776 
8 56 
8 55 
8 53 
8 52 
89 
86 



6 
5 



1303 
1614 
17 21 



4 1304 

5 1664 

6 1622 



1579 
1536 
1545 
14 55 
80 02 
157 4 
1582 
1689 
1267 



3 
4 
1 
1 
1 
1 
6 



5 
6 
6 
5 
5 
6 




44 
30 
15 
16 



83 1893 01 



10 
10 
60 

60 
05 



6 
6 
9 




3 
3 
8 
2 
8 

6 



69 
24 
69 

69 
24 
69 



2 
6 
3 
1 
6 
2 
3 
2 



800 3 
0002 
002 
180 8 
1272 
1832 
1272 
1832 
8 02 

185 4 
001 
1896 
189 6 



1274 
1274 

2 718 
6530 

1872 
5 000 
204 9 




4 
1 
2 
5 



258 
29 4 
818 
222 
714 



6 66 6 



1256 
1267 
1526 

1707 
1267 
1576 



1 

8 


1 

[ 
1 



1 

[ 



28 
00 
00 
56 
B 
27 
00 
25 



3 
2 
3 

1 



183 
184 
1750 
1811 

187 8 

188 7 
187 
188 
144 



4 L 
3>C» 



'L+Mbj for LOG E 
mL+b, for LOG 10 



Normalize and round. 



180 4' 

1 8 6 3 l In case a carry occurs 
1 8 4 f in rounding, correct.* 
180 4. 



183 5 



182 9 
188 1 

183 1 
2818 
659 7 

1874 

908 6 

005 2/ 



509 3 
4819 
1818 
222 2 
285 7 

666 6 
000 1 



1614 
172 1 
1579 

1664 
1622 
1579 



Stop on loss of two digits, 
return to log. progr. on start. 



• Card out of place. 



Constants 



153 

154 

14 

80 

15 

158 

168 

126 



6 

5 

5^ 

2 

4 

2 

9 

7 



C ] 



19. . SINE AND COSINE . 
02-303 . SIN R 



Og-30/i . COS R 



Common 3tepa . 

(this step needed in degree prog.) 



Get B-Bi ,b. Split 
up and store. 











-xviii- 




19 


16 


1575 


68 


8 002 


1538] Test3 in sine routine: 


19 


17 


1538 


15 


1691 


129 5 


19 


18 


1295 


45 


129 8 


17 4 9V Test 47-b. 


19 


1 9 


1298 


46 


156 9 


167 4 J 


19 


30 


1749 


67 


1252 


17 8] 


19 


21 


1708 


15 


1711 


1 6 6 6^ if b-47, test |b|-.0025 


1 9 


3 2 


1666 


44 


156 9 


1674J ' ' 


19 


3 3 


1674 


60 


127 8 


14 4 5 If B < .0025, go to store C« 


19 


24 


1569 


68 


1272 


162 7) ' 


19 


25 


1627 


15 


144 4 


1 6 1 1 If |B| > O 0025, test 49-b, 


19 


2 6 


1601 


3 5 


01 


1 7 9 j go to special routines if 


19 


27 


1709 


44 


1276 


17 1 4J |49- D | > 10. 


19 


28 


1714 


30 


00 5 


17 2 7) 




19 


2 9 


1727 


46 


1682 


1282 




19 


30 


1 2 82 


69 


163 5 


168 81 


For |49-b| < 10, set 


19 


31 


1688 


22 


129 3 


171 5 


shift instruction. 


19 


32 


1682 


69 


153 5 


173 8 




19 


3 3 


1738 


22 


129 3 


174 6. 




19 


3 4 


17 46 


15 


1702 


1 7 1 0) If b > 52, stop to indicate 


19 


3 5 


1710 


46 


137 5 


17 15V loss of accuracy; continue 


19 


3 6 


1 375 


01 




17.1 5 J on start 


19 


37 


1715 


60 


1618 


142 6] 


19 


38 


1426 


19 


1252 


1718 


19 


39 


1718 


35 


2 


12 9 3 Compute the fractional 


19 


40 


1293 


r j ri49~bh 


171 6>part F of B/2tt, then 


19 


41 


1716 


30 


o r 


172 3 


G-F if F ^ 0, G=F+1 if F < 0. 


19 


42 


1723 


46 


1476 


174 1 




19 


43 


1476 


10 


188 6 


174 1 




19 


44 


17 41 


61 


8 00 3 


170 4. 




19 


4 5 


1704 


10 


800 3 


17 13] 


19 


46 


1713 


10 


800 3 


1283Va-4Gfor SIN, 1+4G for COS. 


19 


47 


1283 


u 


188 6 


1742J 


19 


48 


1742 


8001 


174 7llf a £i, a f -a 
170 6) 


19 


4 9 


17 47 


46 


170 5 


19 


50 


1705 


10 


8001 


13 6 2*) if 1< a < 2, a'*2- a 
1748J 


19 


51 


1262 


46 


126 6 


19 


52 


1 366 


10 


8001 


12 7 lllf 2< a<3, a*-2-a 
163 4J 


19 


53 


1271 


46 


1676 


19 


54 


1676 


10 


8 001 


1540 If3<a<5, a'-ar4. 


19 


5 5 


1706 


11 


8 001 


15 4 0) 


19 


56 


1540 


61 


800 3 


17 4 8 [.These angle reductions 


19 


57 


1748 


30 


000 9 


17 2 [carried out. 


19 


58 


1634 


11 


8 001 


174 8 J 


19 


5 9 


1720 


44 


138 3 


17 2 4ilf <x*-l, set S-l. 
12 8 8j 


19 


60 


1724 


20 


128 5 


19 


61 


12 88 


60 


8001 


17 4 4] 
174 5 


19 


62 


1744 


19 


8 001 


19 


6 3 


17 45 


21 


125 


1712 




19 


6 4 


1712 


60 


800 3 


1668 




19 


6 5 


1668 


19 


157 2 


1730 




19 


66 


1730 


60 


800 3 


1387 




19 


6 7 


1387 


11 


1490 


129 6 




19 


6 8 


1296 


19 


1250 


136 8 


>S'0 ] a ,+ 03a ,3 +'»»*0 g a % * 


19 


6 9 


1368 


60 


8 00 3 


172 8 




19 


70 


1728 


10 


153 4 


129 




19 


71 


1290 


19 


1250 


1371 




19 


7 2 


1371 


60 


8 00 3 


128 




19 


7 3 


1280 


11 


168 4 


1389 





-xLx- 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



1 
1 
1 
1 
1 

1 
1 
1 
1 



7 
7 
7 
7 
7 
7 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



80 
81 
82 
83 
8 4 
8 5 
86 
8 7 
8 8 
89 



9 
9 
9 
9 

9 
9 
9 
9 
9 
9 



00 
01 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



9 
5 
4 
4 
6 

2 
9 
3 
2 
4 
5 
9 
7 
3 
5 
7 
6 
6 




256 

526 

7 
6 
9 
1 



7 
6 
1 
1 

1 



35 
35 



19 112 13 53 



1 
3 
6 
1 
1 
6 
3 
4 
3 
3 
3 
3 
6 
3 
1 



9 



9 
5 
5 
4 
1 
5 
1 
5 

6 
1 



46 
10 
11 
46 
60 



2 
6 
6 
1 
6 
4 
9 
3 
3 



15 



07 
64 
15 



1250 
001 
800 3 
139 8 
128 5 
8 00 
000 

38 

00 

00 

00 

00 

00 



3 
1 
3 
3 
2 
2 
1 
2 



8 00 
128 
1 69 
169 
139 
186 



1272 
188 6 
8 02 
188 6 



7 500 




9 

4 
9 
5 
7 



3 
4 
1 
3 
9 
3 
9 



49 

69 138 4 



1 
1 
1 
1 
1 

1 
1 
1 
1 



1739 
150 5 
1585 
1698 
15 55 

13 84 
1385 
1261 
1167 



24 
69 
2 4 
20 
60 

21 
21 
35 
17 



129 3 
126 1 
153 9 
128 9 
116 7 



0002 
4 53 2 



127 
173 
139 
172 
164 
167 
127 
128 
172 
137 
139 
37 



5 

4 

4 

6 



2 

9^ 

4 

2 

9 

5 

9 



Round and normalize, 
1 7 3 7 7" go to store 
12 6 3 



5 
7 
4 4*5 



72 
28 



144 
138 



5J 




Go to stop if b ^ 59, 
set a=0 if b £ 39. 



128 3 

157 5 
174 2 
162 7 
174,2 
157 4 



171 6 

17 16?- Constants 



943 
484 
765 
679 
711 



6318 
005 



173 9 Og-353 . SIN D 



19 113 1354 69 1385 1739 




02-354. COSD 

Common steps 

Set an instr. at the end of 

MFT to get out to SIN R or COS R. 

Go to MET with 6 in place of A 

and 2rr/360 in place of B. 



Constants 



20 



20. ARC TANGENT . 
1 1305 69 1758 1761 2 «-305 . ARC TAN R 



-XX- 



20 



1355 69 1908 1761 0?-355 . ARC TAN D 



Common steps 



20 


3 


1761 


24 


126 4 


1917 


20 


4 


1917 


30 


00 3 


192 5) 


20 


5 


1925 


1 5 


192 8 


800 2 


20 


6 


8002 


60 


129 7 


1941V 


20 


7 


1941 


21 


190 


20 


8 


1900 


3 


002 


175 7 


20 


9 


17 57 


21 


126 5 


176 8J 


20 


10 


1768 


46 


1771 


177 2) 


20 


1 1 


1771 


69 


177 4 


177 8y 


20 


1 2 


1772 


69 


177 5 


177 8 


20 


1 3 


17 78 


24 


128 3 


17 8 6J 


20 


1 4 


1786 


68 


8 002 


179 5) 


20 


1 5 


17 95 


1 5 


156 8 


1773 


20 


16 


1773 


35 


001 


1779 ' 


20 


17 


1779 


44 


178 3 


17 8 4 


20 


1 8 


1783 


4 6 


1937 


193 0. 


20 


19 


1784 


30 


00 5 


189 7) 


20 


20 


1897 


69 


1901 


175 4V 


20 


21 


17 54 


22 


1260 


1763J 


20 


22 


1763 


46 


176 6 


1918 


20 


2 3 


1918 


16 


1921 


192 6) 


20 


2 4 


1926 


46 


192 9 


193 0' 


20 


2 5 


1930 


60 


129 7 


12 6 4J 


20 


2 6 


1929 


65 


126 5 


17 6 91 
1777J 


20 


2 7 


17 69 


69 


1922 


20 


28 


17 66 


60 


1919 


192 3 


20 


2 9 


1923 


64 


126 5 


177 0- 


20 


30 


1770 


69 


1924 


177 7. 


20 


31 


1777 


24 


128 5 


17 8 8' 


20 


3 2 


17 88 


35 


001 


1260. 


20 


3 3 


1260 


30 


{49-b] 
8 2 


17 8 1. 


20 


3 4 


17 81 


67 


17 8 9) 


20 


3 5 


17 89 


20 


129 3 


179 6 


20 


36 


17 96 


1 5 


179 9 


1752) 


20 


3 7 


1752 


4 4 


175 5 


175 6) 


20 


3 8 


17 55 


60 


175 9 


1913) 


20 


3 9 


1913 


19 


129 3 


192 7 


20 


40 


19 27 


3 5 


01 


193 3 


20 


41 


19 33 


10 


166 5 


192 


20 


42 


1920 


21 


127 7 


17 8 


20 


4 3 


1780 


60 


129 3 


179 7 


20 


44 


17 97 


11 


175 9 


176 4 


20 


4 5 


1764 


64 


127 7 


190 2 


20 


46 


19 02 


69 


190 5 


190 9, 


20 


47 


17 56 


65 


129 3 


11 4 71 


20 


48 


1147 


3 5 


01 


17 5 3 1 


20 


49 


17 53 


69 


190 6 


190 9 


20 


50 


1909 


24 


1267 


157 0) 


20 


51 


1570 


20 


127 8 


19 3 1") 


20 


52 


1931 


60 


8 2 


193 9 


20 


53 


19 39 


1 9 


8 1 


194 5 


20 


5 4 


1945 


21 


125 


190 3 


20 


5 5 


19 03 


61 


8 00 3 


1911 


20 


56 


1911 


19 


177 6 


194 


20 


57 


19 40 


60 


800 3 


114 9 


20 


5 8 


1149 


10 


18 5 6 


1762 



Get B-B] ,bj store Bi 



Prepare to use the formula 
arctan (-x) ■ -arctan x 
at the end if B < 0. 



Test 49-b; go to special 
routines if |49-b| > 10 



Prepare for shift of 49-b. 

If b > 49, go to arctan l/x. 
If b ^ 46 » set arctan x ■ x, 
go to end. 

If 46 < b £ 49, go to 

series evaluation 

If b > 49, calculate l/Bj . 



Shift to get x-B or x-l/B 
into fixed decimal form 

If jxl > .28, go to calc. 



z = 



1 + 



z2_. 



x y 



j y=.6 



1- |x|y* 



z « jxj. Prepare to 
add or arctan y. 



-xxi- 



20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 

20 
20 
20 
20 
20 
20 
20 



20 
20 



2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



5 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 



80 
81 



8 
8 
8 
8 
8 
8 



88 

8 9 

9 



1 
2 
3 
4 
5 
6 
7 
8 
9 



00 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



762 
7 9 
9 

1 
3 
4 
6 



8 

4 
8 
8 




8 7 
46 



91 

7 8 

7 

9 

2 

2 

9 

9 



7 
4 
2 
2 
91 
32 



6 

8 

4 

1 

934 
9 4-2 

1 5 

9 

1 

3 

4 

6 

6 

8 

6 

3 

3 



3 
6 
5 
4 
5 
4 
3 
4 
6 
7 



58 
08 
28 
74 



7 

6 



2 

2 

1 

2 

9 

5 

6 





8 

31 

80 

01 

18 

9 4 



1 
6 
1 
1 
6 
1 
1 
6 
1 
1 
3 
6 
1 
1 
6 
1 

[ 
1 

3 

4 

3 

3 

3 

3 



9 

1 
9 


9 

1 
9 



9 



] 
5 

5 

4 

1 

5 

1 

5 



60 
36 
11 
10 

[ ] 



o 1 



65 



21 
60 
61 
60 
49 
30 

65 



6 
9 

1 
1 
1 
2 




125 
8 00 3 
185 5 
1250 
8 00 3 
185 



1 
8 
1 
1 

8 
1 
1 
8 

C 
8 

1 



1 











80 



3 

3 
2 

1 
3 
5 
8 
3 

] 
3 

1 

1 

2 

3 

2 

2 

1 

2 



8 002 
1518 
800 3 

179 4 
1701 



128 3 

800 3 
800 3 



000 3 
8 00 3 
100 
800 3 
2000 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

C 

l 

l 
l 
i 
l 
l 
l 
l 
l 
i 
l 
l 



arc tan z ■ 

n z 2 z 4 z 10 t 

L 3 5 11 J 

(coefficients from LOG) 




(correction at the end) 

arctan |x| -arctan z+£0 or arctan y] 

If b>49, arctan |x| -nr/2- arctan l/|x|. 

If arctan Jb| > 1, round 
in 7th decimal place, if 
< 1, round in 8th. 

Normalize arc tan |b| 



C » arc tan B in radians 

For degrees, go to MPY 

by 360/2tr. 

For b > 59, arc tan |b| - n/2. 



144 5 

193 6 

194 1 
1264 
1264 

178 1 

1910 



1 9 4 3 > Constants 



1731 1285 
1580 128 5 

4041 9 500 



801 



15 7079 6327 
005 
57 2957 8051 
50 



20 114 1143 19 1801 1267 Correction in arc tan z. 
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